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PREFACE 

This  Memorandum  describes  a  model  for  targeting  strikes  against 
a  lines-of-communication  (LOC)  network  and  evaluating  the  effect  of 
carrying  out  the  strategy  that  results.  It  does  not  handle  vehicle 
interdiction,  which  is  treated  in  [4].  The  model  was  developed  while 
the  authors  were  members  of  a  special  RAND-AFGOA  (Operations  Analysis 
Office,  Hq  USAF)  joint  Southeast  Asia  study  group.  Both  RAND  and 
AFGOA  are  currently  applying  the  model  and  will  report  their  results 
to  the  Air  Force  as  part  of  a  follow-on  effort  of  the  group's  activi¬ 
ties.  The  model  generalizes  an  earlier  one  developed  by  Durbin  [l], 
based  on  the  mathematics  of  Wollmer  [6  J  and  applied  in  RAND  and  AFGOA 
studies.  The  mathematical  basis  for  the  model  of  this  Memorandum  may 
be  found  in  [7 ]. 

While  the  model  was  originally  intended  for  immediate  use  by  the 
RAND-AFGOA  study  group  and  later  internal  use  by  both  RAND  and  AFGOA, 
several  other  agencies  have  expressed  interest  in  it  as  an  aid  in  their 
own  research.  Among  these  are  the  Assistant  Chief  of  Staff,  Studies 
and  Analysis,  Hq  USAF;  Army-Air  Force  Intratheater  Transportation  Re¬ 
quirements  Study  Group;  Directorate  of  Force  Planning  Analysis,  Office 
of  the  Chief  of  Staff,  U.S.  Army;  and  the  Boeing  Company.  It  should 
also  prove  useful  to  other  groups  concerned  with  targeting  LOCs,  eval¬ 
uating  LOC  interdiction  effectiveness,  predicting  the  value  of  LOC 
interdiction  tactics,  or  with  constructing  and  repairing  transportation 
systems. 

M.  J.  Ondrasek  is  a  Captain  in  the  Directorate  of  Operations, 

Hq  USAF  (AFXOXSB) . 
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SUMMARY 


This  Memorandum  presents  a  model  for  developing  and  evaluating  a 
targeting  strategy  against  an  opposing  force's  lines  of  communications 
(LOCs)  .  There  are  several  ways  to  measure  how  effective  a  strategy  is 
in  reducing  the  usefulness  of  the  enemy's  LOCs.  Among  these  are  reduc¬ 
tion  in  throughput  (flow),  increase  in  cost  of  meeting  throughput 
requirements,  and  combinations  of  the  two. 

LOCs  are  represented  by  a  network  of  directed  arcs  (line  segments) 
and  nodes  (points).  Arcs  may  represent  road,  rail,  or  waterway  seg¬ 
ments,  or  points  of  transshipment.  They  are  characterized  by  their 
beginning  and  ending  nodes,  upper  and  lower  flow  bounds,  interdicted 
and  uninterdicted  unit  flow  costs,  repair  times,  repair  costs,  and 
probabilities  that  attempted  strikes  result  in  successful  interdiction. 
Nodes  are  junctions  where  two  or  more  arcs  intersect. 

The  model  operates  in  daily  cycles,  the  user  specifying  the  number 
of  days  and  strikes  per  day.  At  the  start  of  each  day,  certain  inter¬ 
dicted  arcs  are  returned  to  operation  (as  specified  by  their  repair 
times  and  last  time  struck);  the  capacities  and  unit  flow  costs  of 
these  arcs  are  restored  to  uninterdicted  values.  Then  the  appropriate 
number  of  strikes  are  targeted,  one  at  a  time,  against  the  LOC  arcs  by 
the  second  algorithm  in  Ref.  [7].  The  algorithm  chooses  targets  based 
on  the  cost  and  time  involved  in  repairing  the  arc  and  on  how  the  arc's 
disablement  will  reduce  network  effectiveness.  The  procedure  is  optimal 
when  one  strike  is  targeted,  and  approximates  optimality  when  more  than 
one  strike  is  targeted.  At  the  end  of  each  strike,  total  throughput  and 
the  necessary  cost  to  the  LOC  user  of  achieving  this  throughput  are  pro¬ 
duced  as  output. 

Two  other  outputs  are  available  to  the  model  user,  and  either  or 
both  may  be  requested  after  each  strike  or  at  the  end  of  the  day.  The 
first  is  a  detailed  status  of  the  entire  network,  including  individual 
arc  flows;  the  second  is  a  profile  of  total  flow  versus  required  user 


cost . 
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I.  INTRODUCTION 


This  Memorandum  presents  a  model  for  targeting  strikes  against  an 
opponent's  lines  of  communication  (LOCs)  and  evaluating  their  effective¬ 
ness  in  reducing  his  ability  to  deliver  supplies.  Effectiveness  of 
strikes  may  be  measured  by  reduction  in  throughput,  cost  of  operation 
to  the  LOC  user,  or  combinations  of  these.  Operating  costs  may  be 
dollar  costs,  manpower  or  personnel  costs,  or  other  appropriate  measures. 

Single-mode  LOCs  or  multimode  LOCs  with  transshipment  between  modes 
can  be  represented  by  a  network  of  arcs  and  nodes.  Arcs  may  represent 
road,  rail,  or  waterway  segments,  or  points  of  transshipment  between 
modes.  Nodes  are  junctions  where  two  or  more  arcs  intersect.  Arcs  are 
characterized  by  their  beginning  and  ending  nodes,  lower  and  upper 
bounds  on  flow,  and  unit  flow  costs.  In  this  model,  an  arc's  upper 
bound  on  flow  and  unit  flow  cost  are  parameters  that  are  functions  of 
the  number  of  strikes  targeted  against  the  arc.  These  functions  de¬ 
pend  for  their  arguments  on  the  values  of  these  parameters  when  the 
arc  is  struck  successfully,  when  unstruck  or  struck  unsuccessfully, 
and  on  the  probability  it  has  been  struck  successfully;  it  is  the  latter 
quantity  that  varies  with  the  number  of  strikes. 

Many  network  flow  problems  can  be  formulated  to  find  a  minimum 

* 

cost  circulation  flow  by  artificial  devices.  Among  these  are  maximiz¬ 
ing  flow  from  a  source  node  to  a  sink  node,  meeting  a  required  flow 
between  two  nodes  at  minimum  operating  cost,  and  combinations  of  these 
two.  It  is  assumed  that  the  effectiveness  of  the  LOCs  to  the  user  may 
be  formulated  in  this  context.  The  effectiveness  of  a  strike  is  mea¬ 
sured  in  terms  of  how  it  reduces  LOC  usefulness  to  its  user. 

The  programed  model  operates  in  daily  cycles.  Each  day  a  given 
number  of  strikes  are  targeted  against  the  LOC  arcs  after  certain  arcs 
from  among  those  previously  struck  are  returned  to  their  uninterdicted 


As  will  be  shown  later,  this  cost,  obtained  partly  by  artificial 
devices,  often  differs  from  actual  operating  costs. 
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condition.  Arcs  are  selected  as  targets  on  the  basis  of  what  strikes 
on  them  contribute  towards  reducing  LOC  effectiveness  to  its  user. 
Specifically,  the  arc  chosen  for  the  next  strike  is  the  one  that  will 
maximize  the  repair  cost  plus  the  product  of  the  repair  time  and  the 
cost  increase  of  a  minimum  cost  circulation  flow. 
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II.  NETWORK  REPRESENTATION  OF  LOCs 


SINGLE  MODE 

A  single-mode  system  of  transportation  routes  such  as  roads  or 
waterways  can  be  represented  by  a  network  of  nodes  and  directed  arcs. 
Specifically,  for  a  highway  system,  a  node  is  the  intersection  of  two 
or  more  roads,  and  an  arc  represents  a  road  segment  joining  two  nodes. 

For  a  one-way  segment,  the  beginning  node  is  the  point  where  the  traf¬ 
fic  enters  the  segment  and  the  ending  node  the  point  where  traffic  leaves 
it.  A  two-way  road  segment  can  be  represented  by  two  arcs,  one  in  each 
direction.  Similar  representations  hold  for  railroads,  waterways,  and 
other  modes.  As  an  example,  consider  the  network  below. 


a 


Fig.  1  —  A  simple  network 

Figure  1  might  represent  a  road  system  with  segments  joining  locations 
S  and  a,  a  and  Tt  b  and  T,  a  and  T,  d  and  Tt  S  and  a3  S  and  dt  a  and  b, 
b  and  c}  and  a  and  d.  All  roads  are  two-way,  and  hence  each  segment  Is 
represented  by  two  arcs. 

Each  arc  has  certain  parameters:  flow  capacity  (the  upper  bound 
on  flow),  lower  bound  on  flow,  and  cost  per  unit  flow.  As  an  example, 
let  us  suppose  that  the  parameters  for  the  arc  in  Fig.  1  joining  r>  and 


a  and  In  die  direction  of  to  a,  designated  by  the  symbol  are 

as  in  Table  1. 


Table  1 

PARAMETERS  FOR  ARC  ( S,a )  IN  FIG.  1 


Parameter  Value 

Flow  .  6 

Upper  Bound  .  8 

Lower  Bound  .  0 

Coat  .  3 


In  Table  1,  the  actual  flow  from  S  to  a  is  6.  This  can  be  tons/ 
day,  trucks/hour  or  any  other  appropriate  measure.  Other  values  of  flow 
are  possible,  but  only  between  the  inclusive  llmitB  of  aero  and  8.  Lower 
bounds  normally  are  zero,  as  in  this  example;  but  as  we  shall  see,  there 
is  sometimes  a  reason  for  having  them  greater  than  zero.  The  uppet  bound 
here  could  represent  the  road's  physical  capacity  to  handle  traffic.  The 
cost  per  unit  flow  of  3  is  in  whatever  units  are  used  for  the  example; 
possibilities  include  dollars  per  ton  moved  from  S  to  a,  man-hours  per 
vehicle  traveling  from  S  to  a,  and  vehicle  hours.  The  total  cost  for 
the  flow  of  6  on  this  arc  would  be  3  times  6,  or  18.  Note  that  one  would 
normally  expect  that  except  for  flow,  all  parameters  for  the  arc  in  the 
opposite  direction  (a,5)  would  be  identical  to  those  of  (Sta). 

MULTIMODE  WITH  TRANSSHIPMENT 

A  system  with  several  modes  of  transportation  can  be  represented 
by  first  constructing  a  separate  network  for  each  mode  and  then  connect¬ 
ing  these  networks  at  transshipment  points.  For  example,  lei  us  take 
the  simple  railroad  network  of  Fig.  2,  together  with  the  road  network 
of  Fig.  1.  The  points  S,Tto,  and  d  in  Fig.  2  represent  the  same  points 
os  these  some  symbols  in  Fig.  1.  Suppose  transshipment  is  allowed  only 
at  nodes  5  and  d.  Then  the  combined  road-rail  network  may  be  represented 
as  in  Fig.  3.  Note  that  o  and  o'  (the  two  nodes  representing  point  o ) 
and  T  and  T'  are  not  connected.  This  is  because  flow  entering  point  o 


Fig.  2  —  Network  representation  of  railroad  routes 


a 


d' 

Fig.  3  —  Combined  road-rail  network  of  Fig.  1  and  Fig.  2 


or  point  T  by  one  mode  cannot  leave  by  another.  This  also  means  that 
a  and  c '  (and  similarly  T  and  T')  cannot  be  represented  by  one  node,  if 
the  network  is  to  have  the  property  that  flow  entering  a  node  via  a  par¬ 
ticular  arc  can  leave  that  node  via  any  arc  emanating  from  that  node. 

The  flow,  upper  bound,  lower  bound,  and  cost  on  arc  (d,d')  represent 
the  amount  transshipped  from  road  to  rail  at  that  point,  capacity  of 
the  transshipment  facility,  lower  bound  on  use  of  facility  (usually 


zero),  and  cost  per  unit  of  flow  transshipped.  Flow  on  ( d' ,d )  repre¬ 
sents  the  amount  transshipped  from  rail  to  road,  and  other  parameters 
have  the  same  meaning  as  in  arc  (dtd') . 

SOURCES.  SINKS.  AND  FLOWS 

Very  often  one  is  interested  in  sending  units  of  flow  from  one 
node  to  another.  As  an  example,  one  may  wish  to  transport  equipment 
from  a  supply  depot  to  a  battlefront.  In  such  cases  the  node  the  flow 
originates  from  is  called  the  source  and  the  node  it  terminates  at 
is  called  the  sink. 

As  an  example,  suppose  in  Fig.  1  the  user  wished  to  send  units  of 
flow  from  S  to  T.  Then  5  would  be  the  source  and  T  would  be  the  sink. 

If  only  one  unit  were  sent  from  S  to  T,  a  possible  route  is  the  path 
Sta,b,Ti  this  would  be  represented  by  a  unit  of  flow  on  arcs  ( 5, a )■ 

( a,b ),  and  (btT) ,  and  zero  flow  on  all  other  arcs.  A  flow  of  one  unit 
over  an  alternate  route,  StctT  would  be  represented  by  one  unit  on  (S,c), 
one  on  (c,T) ,  and  zero  units  on  all  other  arcs.  Note  that  in  each  case, 
the  net  flow  out  of  S  and  the  net  flow  into  T  both  equal  one,  while  the 
net  flow  into  or  out  of  any  other  node  is  zero.  In  general  networks, 
if  x  units  were  to  be  sent  from  a  source  node  to  a  sink  node,  then  any 
flow  pattern  accomplishing  this  would  have  net  flows  of  x  units  out  of 
the  source,  x  units  into  the  sink,  and  zero  into  (or  out  of)  any  other 
node.  The  reverse  also  follows;  specifically,  any  flow  pattern  having 
net  flows  of  x  units  out  of  the  source,  x  units  into  the  sink,  and  zero 
units  of  net  flow  into  or  out  of  any  other  node,  can  be  broken  up  into 
units  of  flow  on  routes  from  source  to  sink  totaling  x.  In  addition  to 
these  routes,  however,  there  may  be  cycles  or  closed  loop  routes. 

In  actual  flow  situations,  there  may  be  several  sources  and  several 
sinks.  Such  problems  may  be  handled  in  the  LOC  model  provided  any  source 
■ay  supply  any  sink.  Of  course,  this  condition  is  always  met  if  there 
are  several  sources  and  one  sink  or  one  source  and  several  sinks.  The 
reason  for  this  condition  is  that  it  allows  one  to  convert,  through  arti¬ 
ficial  devices,  the  multisource,  multisink  problem  to  a  single-source, 
single-sink  problem,  and  technically  we  require  but  one  source  and  one 
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slnk.  The  way  this  conversion  is  accomplished  is  to  connect,  with  arti¬ 
ficial  arcs,  all  sources  to  an  artificial  super  source  and  all  sinks  to 
an  artificial  super  sink.  For  example,  suppose  in  the  network  of  Fig.  1 
that  a  and  .  '  were  sources  and  < 1  and  7  were  sinks.  Then  this  network 
could  be  represented  by  the  single-source,  single-sink  network  of  Fig.  4. 


Fig.  4  —  Equivalent  network  to  double  source  and  double 
sink  network  of  Fig.  1.  The  arrows  indicate 
the  direction  of  the  arc.  Lines  without  arrows 
indicate  arcs  in  both  directions. 


The  artificial  or  super  source,  3  (in  Fig.  4),  would  be  the  source 
and  T  the  sink.  Flow  on  the  artificial  arc  (3, a)  equals  the  rate  that 
goods  leave  the  true  source,  a.  The  lower  bound  on  (3, a)  would  normally 
be  zero;  the  upper  bound  might  be  essentially  infinite  or,  if  it  repre¬ 
sented  a  facility  such  as  a  warehouse,  could  be  the  maximum  supply 
rate  at  that  facility;  and  cost,  if  not  zero,  could  be  equal  to  the  unit 
cost  of  using  facilities  represented  by  a •  The  same  interpretation  may 
be  given  to  the  parameters  of  An  artificial  arc  connecting  a 

sink  to  the  artificial  super  sink  node  carries  a  flow  equal  to  the  quan¬ 
tity  per  unit  time  entering  the  sink  it  represents,  a  lower  bound  equal 


r 
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to  zero  or  the  sink  demand,  an  upper  bound  essentially  infinite  or  equal 
to  sink  demand,  and  a  cost  either  zero  or  equal  to  the  unit  cost  of  using 
the  sink  facility.  Note  that  such  a  formulation  as  in  Fig.  4  assumes 
any  Bource  can  supply  any  sink  in  any  proportion.  Thus  this  formulation 
would  not  be  valid  if  it  were  required  that  flow  emanating  from  source  a 
must  terminate  at  sink  d. 
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III.  USER  EFFECTIVENESS  AND  CIRCULATION  FLOWS 

TRANSPORT  USER  EFFECTIVENESS  OF  LOCs 

As  will  be  shown  in  Che  next  section,  strike  effectiveness  is 
based  quite  heavily  on  the  degree  to  which  it  reduces  the  usefulness 
of  the  LOCs  to  the  user.  User  effectiveness  of  LOCs  may  he  measured 
in  several  ways.  Among  these  are  the  following: 

1.  Maximum  flow  from  one  node  (source)  to  another  (sink). 

2.  Minimum  cost  of  meeting  a  required  source-sink  flow. 

3.  Minimum  cost  of  required  flow  if  required  flow  can  be  met; 
otherwise,  the  least  cost  of  maximizing  flow. 

The  first  of  these  is  realistic  when  the  user  is  physically  limited 
by  the  ability  of  his  transport  system  to  handle  traffic.  The  second 
is  reasonable  when  he  is  not  limited  by  his  transport  system  but  is 
constrained  by  a  resource  such  as  dollars  or  man-hours,  or  merely  de¬ 
sires  to  use  resources  efficiently.  The  third  is  one  of  several  pos¬ 
sible  combinations  of  the  first  two,  and  applies  when  the  user  does 
not  know  whether  or  not  he  is  constrained  by  his  transport  system  and 
also  desires  to  use  resources  efficiently. 

CIRCULATION  FLOWS 

A  minimum-cost  circulation  flow  pattern  in  a  network  can  be  found 
with  the  "Out-of-Kilter"  algorithm  Ll].  Such  a  network  problem  requires 
minimizing  total  cost  subject  to  the  constraint  that  all  arc  flows  are 
between  their  upper  and  lower  bounds,  and  the  flow  entering  a  particular 
node  must  equal  the  flow  leaving  it.  Not  so  obvious  is  the  fact  that 
each  user  effectiveness  criterion  discussed  above  can  be  formulated  as 
a  minimum-cost  circulation  flow  problem. 

Converting  source-sink  flows  to  circulation  flows  is  accomplished 
as  follows.  Connect  a  universal  arc  directed  from  the  sink  to  the 
source  and  assign  it  a  flow  equal  to  that  of  the  total  source-sink 
flow.  Letting  5  be  the  source  and  T  the  sink  in  Fig.  1,  this  network 
would  be  redrawn  as  in  Fig.  5. 


-10- 


cl 


Fig.  5  —  letwork  of  Fig.  1  with  artificial  arc  added 


Any  flow  pattern  sending  x  total  units  of  flow  from  S  to  T  for  the 
network  of  Fig.  1  can  be  represented  by  a  circulation  flow  pattern  in 
the  Fig.  5  network  by  assigning  a  flow  of  x  units  to  arc  ( T,S )  and  flows 
on  other  arcs  identical  to  those  of  the  Fig.  1  network.  The  reverse  is 
also  true. 

Network  effectiveness  to  the  LOC  user,  as  measured  by  any  of  the 
three  criteria  of  the  last  subsection,  may  be  found  by  solving,  respec¬ 
tively,  the  following  three  problems. 

1.  Maximize  flow  from  one  node  (source)  to  another  (sink), 

2.  Meet  a  given  or  required  source-sink  flow  at  minimum  cost, 

3.  If  it  is  possible  to  meet  a  required  source-sink  flow,  do  so 
at  minimum  cost.  If  this  is  impossible,  maximize  flow  in  the 
least  costly  way. 

These  problems  may  be  solved  by  suitable  choices  for  parameters 
(other  than  flow)  of  the  universal  arc.  The  problem  of  maximizing  source- 
sink  flow  is  converted  to  one  of  finding  a  minimum  cost  circulation  flow 
by  assigning  a  cost  of  minus  infinity  (or  a  negative  number  whose  abso¬ 
lute  value  is  at  least  as  great  as  the  sum  of  costs  on  all  other  arcs). 
This  insures  that  all  flow  cycles  incorporating  the  artificial  arc  have 
a  negative  cost.  Hence  the  minimum  cost  circulation  problem  is  essen¬ 
tially  one  of  maximizing  flow  on  the  universal  arc,  which  is  the  same 


as  maximizing  source-sink  flow.  The  lower  bound  on  the  universal  arc 
would  he  zero  and  the  upper  bound  infinity. 

A  required  flow,  /•,  at  minimum  cost  may  be  met  by  setting  the  upper 
and  lower  bounds  on  the  universal  arc  equal  to  v.  Cost  on  t lie  universal 
arc,  which  Is  arbitrary,  would  normally  be  set  to  zero. 

The  third  problem  may  be  formulated  by  assigning  the  universal  arc 
a  cost  of  minus  infinity,  a  lower  bound  of  zero,  and  an  upper  bound  of 
t lie  required  flow,  r.  This  assures  the  maximization  of  flow  up  to  a 
value  of  r.  Among  all  flow  patterns  accomplishing  this,  the  least  cost 
one  will  be  selected.  Theoretically,  this  formulation  is  also  appropriate 
when  it  is  known  that  maximum  flow  can  be  reduced  below  required  user 
flow.  For  this  case.  .towever,  the  strike  locations  found  by  the  model 
will  normally  be  closer  to  optimal  if  formulated  in  the  context  of  the 
first  problem. 


IV.  MODEL  OPERATION 


DESCRIPTION 

From  a  mathematical  viewpoint,  it  is  assume’  that  the  effectiveness 
of  the  LOC  network  to  its  user  may  be  measured  in  terms  of  how  cheaply 
he  may  achieve  a  circulation  flow.  This,  of  course,  includes  the  phys¬ 
ical  effectiveness  criteria  discussed  in  Sec.  Ill,  namely,  maximum  flow, 
minimum  cost  of  required  flow,  and  the  combined  criterion.  The  model 
attempts  to  reduce  LOC  effectiveness  over  time,  taking  into  account  re¬ 
pair  cost,  and  consequently  tries  to  maximize  cost  by  targeting  strikes 
against  the  arcs  of  the  LOC  network. 

Striking  an  arc  decreases  its  capacity  and  increases  its  unit  flow 
cost  for  a  specified  period  of  time,  after  which  its  parameters  may  be 
restored  to  their  original  value  by  incurring  a  repair  cost.  This  spe¬ 
cified  period  is  referred  to  as  the  repair  time.  In  targeting  a  strike, 
the  arc  selected  will  be  the  one  that  maximizes  repair  cost  plus  the 
product  of  repair  time  and  cost  increase  of  a  minimum  cost  circulation 
flow.  Specifically,  if  k  is  the  cost  of  a  minimum  cost  circulation  flow 
before  allocating  a  strike,  k..\  this  cost  after  allocating  a  strike  to 

arc  t.  .  and  r..,  respectively,  the  repair  time  and  repair  cost 

I'J  l  J 

for  arc  (t,j)  ;  then  the  strike  will  be  targeted  against  the  arc  that 
maxi mizes 


(1) 


(k..-k)t.  . 
ij  iJ 


+  r. 


i«7 


The  above  function  is  quite  general  and  one  may  derive  special  cases. 
For  example,  repair  costs  may  be  eliminated  from  the  calculation  by  set¬ 
ting  them  equal  to  zero.  The  immediate  effect  of  striking  arcs  (l.e., 
not  taking  time  into  consideration)  may  be  maximized  by  setting  all  re¬ 
pair  times  equal  to  one  and  letting  the  campaign  run  for  one  day. 


ARC  PARAMETERS  AND  STRIKES 

All  strikes  allocated  against  a  given  arc  are  directed  against  a 
* 

common  target.  Each  strike  has  identical  but  independent  probabilities 
* 

To  see  how  this  restriction  may  be  removed,  see  Ref.  7. 
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of  success.  True  arc  capacities  and  unit  flow  costs  are  double  valued. 

The  first  of  these  values  for  each  are  uninterdicted  values  that  apply 
to  an  arc  when  either  no  strikes  are  directed  against  it  or  some  strikes 
are  directed  against  it  but  all  are  unsuccessful;  the  second  are  inter¬ 
dicted  values  that  apply  when  the  arc  has  been  struck  successfully.  Of 
course,  interdicted  capacities  (unit  flow  costs)  are  lower  (higher)  than 
uninterdicted  ones. 

The  model  expresses  arc  capacities,  unit  flow  costs,  and  repair 

times  as  functions  of  the  number  of  strikes,  based  on  expected  values. 

When  possible,  steady  state  movement  is  such  that  all  flow  across  struck 

arcs  will  occur  when  strikes  are  unsuccessful  or  when  the  arc  has  been 

repaired.  Thus  a  targeted  arc  is  either  used  to  its  full  capacity  when 

uninterdicted  or  not  at  all  when  interdicted.  Furthermore,  it  is  assumed 

that  men  or  equipment  delayed  in  transit,  which  figure  in  unit  flow  costs, 

can  always  be  used  productively  for  other  tasks.  Thus  no  cost  penalty  is 

assessed  for  time  spent  waiting  for  arcs  to  be  repaired.  For  arc  ( ), 
INT  UNINT 

let  u . .  and  u..  represent  interdicted  and  uninterdicted  capacities 

INT  UNINT 

or  upper  flow  bounds,  a.,  and  o.  .  interdicted  and  uninterdicted 
INT 

unit  flow  costs,  r..  the  true  interdicted  repair  time,  and  p.  .  the 

g  gJ 

probability  of  a  strike  being  successful;  and  then  let  u .  .,  e  .  and 

8  t<7  X ^ 

v.  .  represent,  as  seen  by  the  model,  upper  flow  bound,  unit  flow  cost, 

and  repair  cost  for  arc  (£,j)  when  s  strikes  are  targeted.  Then  the 
following  relations  hold; 


(2) 


8  ,8,  INT  „  ,8  UNINT 

U..=  fl-(l-p.J  )u..+  (L-p.J  U.  . 

IJ  T'J 


(3) 


8 

a .  . 
T'J 


UNINT  .  ...  _  ,,  .8,  UNINT 

io..  for  flow  units  up  to  (1-p.  J  u.  . 

T'J 


■  INT  .  .  ...  .8  UNINT 

[a..  for  flow  units  beyond  (1-p.J  u.. 

1-J  tj 


8  /I  /I  1 6  i  INT 
r.  .  =  (1- ( 1-p .  J  )r. .  . 

tj  ij 


(4) 
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Repair  costs  in  (4)  are  cumulative.  Thus,  if  arc  (i,J)  already 

has  been  struck  a  times,  t he  p ,  .  used  in  (1)  would  be  the  difference 
i'+l  i‘ 

between  r  and  r.  ..  Repair  time,  t  .  is  the  same  for  all  levels 

i.l  1 

of  strike.  The  quantities  k  and  k..  are  both  costs  of  minimum  cost 
circulation  flows,  but  k  is  the  value  that  results  when  capacity  and 
unit  flow  cost  for  arc  (i,j)  are  u  .  .  and  a,.,  respectively,  while  k . 


*1 


',+ 1 


n+1 


t,l 


results  when  these  quantities  are  u".  .  and  c‘.  .  .  Of  course,  the  ca- 

i.]  i,l 

pacity  and  unit  flow  costs  for  arcs  other  than  ( i,,j )  are  the  same  for 

the  calculation  of  k .  .  as  they  are  for  k  and  are  those  which  correspond 

i-J 

to  the  number  of  strikes  already  targeted  against  them. 


OPERATION 

The  model  operates  in  daily  cycles,  the  model  user  specifying  the 
number  of  days  and  strikes  per  day.  At  the  start  of  each  day,  capacities 
and  unit  flow  costs  on  arcs  scheduled  to  be  returned  to  their  unstruck 
condition  (as  specified  by  their  repair  times  and  last  time  struck)  are 
restored  to  uninterdicted  values.  Then  the  appropriate  number  of  strikes 
are  targeted,  one  at  a  time,  against  the  LOC  arcs  by  the  second  algo¬ 
rithm  in  Ref.  [7].  Each  time  the  arc  chosen  to  be  struck  is  one  for 
which  expression  (1)  is  a  maximum.  Note  that  this  is  a  marginal  allo¬ 
cation  procedure,  thus  giving  solutions  with  the  property  that  if  an 
arc  is  struck  when  x  number  of  strikes  are  targeted,  it  will  also  be 
struck  when  more  than  x  strikes  are  targeted.  Since  optimal  placement 
of  strikes  would  not  necessarily  have  this  property,  the  strike  locations 
are  optimal  for  one  strike  but  only  approximate  optimality  for  multiple 
strikes . 

After  each  strike,  information  on  the  arc  struck,  total  source- 
sink  flow,  and  total  cost  of  this  flow  is  always  output.  Since  the 
cost  on  the  universal  arc  is  not  really  incurred,  but  merely  a  device 
to  convert  a  particular  type  of  source-sink  flow  problem  to  one  of  cir¬ 
culation  flows,  it  is  not  considered  in  the  printed  output  cost  calcu¬ 
lation.  Costs  on  artificial  arcs  are  also  excluded  from  this  calcula¬ 
tion;  however,  the  user  may  elect  to  include  these  costs  by  inputting 
them  as  real  arcs.  In  this  regard,  the  cost  given  as  output  differs 
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from  the  quantity,  k,  referred  to  in  expression  (1).  In  particular, 
targeting  an  additional  strike  may  result  in  the  output  of  a  lower  total 
cost  (and  source-sink  flow),  but  could  never  result  in  a  lower  k  value. 

In  addition,  a  detailed  status  of  the  network  and  arc  flows  is  printed 
at  the  beginning  of  each  day  and  at  the  end  of  the  campaign. 

Two  other  output  features  are  also  available  at  the  user's  option. 
First,  the  status  of  the  entire  network  may  be  printed  after  each  strike. 
Second,  a  profile  of  cost  versus  source-sink  flow  may  be  obtained  either 
at  the  beginning  of  each  day  and  after  the  campaign,  or  at  the  beginning 
of  each  day  and  after  each  strike. 

When  the  second  option  is  required,  total  source-sink  flow  is  di¬ 
vided  by  the  specified  number  of  profile  points.  Multiples  of  this  num¬ 
ber  then  serve  as  profile  points.  The  cost  of  meeting  flows  is  found 
explicitly  at  each  point.  For  values  of  source-sink  flow  lying  between 
profile  points,  cost  may  be  approximated  by  interpolation. 

It  is  important  to  note  that  network  status  and  profile  points  at 
the  beginning  of  a  day  may  differ  from  those  at  the  end  of  the  last 
strike  of  the  previous  day.  The  reason  for  this  is  that  in  the  former 
case  any  arcs  scheduled  to  return  to  the  f  r  unstruck  condition  are  re¬ 
stored  to  their  unstruck  values,  while  in  the  latter  case  tech  resto¬ 
ration  has  not  yet  taken  place. 
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INPUT  DATA 

Control  Data 

1.  N 

2.  PERIOD 

3.  ICLASS 


A.  10P 


5.  I ART 


Appendix  A 

PREPARATION  OF  INPUT  DECK 


is  the  number  of  strikes  to  be  targeted  each  day. 
is  the  number  of  days  for  which  the  model  is  to  target 
strikes  and  evaluate  network  status. 

specifies  the  classification  to  be  printed  on  the  out¬ 
put.  The  code  is 

1.  Unclassified 

2.  Confidential 

3.  Secret 

4.  Top  Secret 

5.  Confidential,  no  foreign  dissemination 

6.  Secret,  no  foreign  dissemination 

7.  Top  Secret,  no  foreign  dissemination 

8.  For  Official  Use  Only 

controls  the  output.  If  IOP  is  1,  a  paragraph  is  written 
on  the  results  of  each  strike,  and  a  detailed  statement 
of  network  status  is  given  at  the  beginning  of  each  day 
and  after  the  campaign.  If  IOP  is  2,  a  paragraph  on  the 
results  of  each  strike  and  a  detailed  status  of  the  net¬ 
work  is  given  after  strike  as  well  as  at  the  beginning 
of  each  day.  If  IOP  is  equal  to  3,  the  information  printed 
is  the  same  as  when  it  is  1,  except  that  the  results  of 
a  day's  strikes  are  written  in  tabular  form  rather  than 
in  a  series  of  paragraphs.  When  the  campaign  begins,  the 
network  status  printout  consists  of  information  for  all 
arcs  so  that  one  has  a  record  of  the  input  network  struc¬ 
ture;  however,  subsequent  listings  contain  only  arcs  with 
non-zero  flow. 

is  the  number  of  artificial  arcs  in  the  network  (i.e., 
arcs  connecting  true  sources  to  the  super  source  and 
true  sinks  to  the  super  sink).  The  universal  arc,  which 
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connccts  t lie  super  sink  to  the  super  source,  is  not  in¬ 
cluded  in  this  count. 

6.  IPROFL  determines  if  and  when  the  profile  of  flow  versus  cost 

is  calculated  and  printed.  If  1,  the  profile  option  is 
used  at  the  end  of  each  day;  if  2,  it  is  used  after  each 
strike;  if  3,  it  is  not  used  at  all. 

7.  NPPTS  is  the  number  of  profile  points  (not  to  exceed  25)  de¬ 

sired.  Of  course,  this  parameter  only  applies  when  IPROFL 
equals  1  or  2 . 

8.  TOWAY  All  arcs  are  considered  one-way  if  this  equals  1.  If 

this  equals  2,  then  all  arcs  with  return  arcs  are  con¬ 
sidered  to  be  two-way  (i.e.,  both  and  (j3i)  exist 

and  are  designated  by  the  same  arc  name  or  as  defined 
below,  ARCNM) .  Whenever  a  two-way  arc  is  struck,  the 
strike  applies  to  the  return  arc  as  well. 

9.  JOBS  is  equal  to  one. 

Arc  Data 

1.  ARCNM  is  the  name  of  the  arc.  For  two-way  arcs,  this  must  be 

the  same  for  both  of  the  directed  arcs.  Otherwise,  it  is 
used  only  for  identification  purposes. 

2.  FROM  is  the  beginning  node  of  the  arc.. 

3.  TO  is  the  ending  node  of  the  arc. 

4.  ICCP  is  the  uninterdicted  capacity. 

5.  LCAP  is  the  interdicted  capacity. 

6.  ICCST  is  the  uninterdicted  unit  flow  cost. 

7.  UCOST  is  the  interdicted  unit  flow  cost. 

8.  REPR  is  the  cost  of  repairing  an  interdicted  or  successfully 

struck  arc. 

9.  ITME  is  the  time  required  for  repairing  an  interdicted  arc. 

10,  LWER  is  the  lower  bound  on  arc  flow. 

11.  P(I, J)  is  the  probability  that  a  strike  is  successful. 
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INPUT  FORMAT 


A 

Denotes 

alphanumeric  data. 

I 

Denotes 

an  integer  right-adjusted  in  its  field. 

D 

Denotes 

a  decimal  fraction.  A  decimal  point  must  be 

included 

in  the 

input . 

First  Card 

Col 

1-78 

Title  of  Run 

Second  Card 

Col 

3-  5 

Number  of  Strikes  (N) 

(I) 

Col 

8-10 

Number  of  Days  (PERIOD) 

(I) 

Col 

13-15 

Classification  Number  (ICLASS) 

(I) 

Col 

18-20 

Print  Mode  (I0P) 

(I) 

Col 

23-25 

Number  of  Artificial  Arcs  (IART) 

(I) 

Col 

28-30 

Profile  Option  (IPROFL) 

(I) 

Col 

33-35 

Number  of  Profile  Points  (NPPTS) 

(1) 

Col 

38-40 

Two-way  Arc  Option  (TOWAY) 

(1) 

Col 

43-45 

Number  of  Jobs  (JOBS) 

(I) 

Third 

Card 

Col 

1-  4 

The  word  "ARCS" 

(I) 

Arc  Data  Cards 


Each  arc  is  defined  by  two  data  cards.  The  first  contains  un¬ 
interdicted  values,  the  second  contains  interdicted  values.  The  cards 
for  the  universal  arc  must  come  first,  those  for  all  artificial  arcs 
next,  and  those  for  all  real  arcs  last. 


First  Arc  Card 


Col 

3-  6 

Name  of  Arc  (ARCNM) 

(A) 

Col 

10-13 

Beginning  Node  (FROM) 

(A) 

Col 

17-20 

Ending  Node  (TO ■ 

(A) 

Col 

22-27 

Uninterdicted  Capacity  (ICCP) 

(1) 

Col 

32-37 

Uninterdicted  Cost  (ICCST) 

(I) 

Col 

42-47 

Repair  Cost  (REPR) 

(I) 

Col 

52-57 

Repair  Time  (ITME) 

(I) 

Col 

62-67 

Lower  Bound  (LWER) 

(I) 

Col 

72-75 

Kill  Probability  P(I,J) 

(D) 
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Second  Arc 

Card 

Col 

3-  6 

Name  of  Arc  (ARCNM) 

(A) 

Col 

10-13 

Beginning  Node  (FROM) 

(A) 

Col 

17-20 

Ending  Node  (TO) 

(A) 

Col 

22-27 

Interdicted  Capacity  (LCAP) 

(I) 

Col 

32-37 

Interdicted  Cost  (UCOST) 

(I) 

Col 

42-47 

Repair  Cost  (REPR) 

(I) 

Col 

52-57 

Repair  Time  (ITME) 

(I) 

Col 

62-67 

Lower  Bound  (LWER) 

(I) 

Col 

72-75 

Kill  Probability  P(I,J) 

(D) 

Next 

Card 

Col 

3-  9 

The  word  "COMPUTE” 
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main 

C  AACNN  .  TMt  NAME  Of  AN  AHC . 

C  Amo .  TMt  SANE  AS  MAXFLOII)  BUT  IN  At  At  NODE. 

C  BSIAAC .  THE  SUBSCRIPT  OF  TMt  Bt ST  AHC  TO  STRlKk. 

C  OAY  . USEO  TO  CUUNI  THE  CAMPAIGN  DAYS. UP  TO  PEPIOO. 

C  OH  TA  .  THE  VALUE  USEO  TO  LONER  THE  UPPER  BIJUNO. 

C  ERROR  .  A  VALUE  SENT  TO  THE  ERROR  SUBNOUTINL  SO  THAT  THE  CORRECI 

C  MESSAGE  IS  PRINTEO. 

C  FLOWIII  . ELOM  ON  ARC  I. 

C  FHUMIII  .  II  AT  INPUT  TIRE  CONTAINS  THE  LOCATION  OF  THE  FREOECESSOR 

C  NODE  IN  GEOREF  OR  UTM. 

C  21  AFTER  NOOES  ARE  NUMBERED  IT  CONTAINS  THE  UNIQUE 

C  NUOE  NUMBER. 

C  I  AH  T  .  THE  NUMBER  OF  ARTIFICIAL  ARCS  IN  THE  NETNORR. 

C  IBSIAC .  IF  THE  BEST  ARC  TO  STRIKE  NAS  TNO-NAY,  IT  CONTAINS  THE 

C  SUBSCRIPT  Of  THE  RETURN  ARC. 

C  ICAPIII  . UNINTEMDIC  TEO  CAPACITY  FOR  ARC  I. 

I  ICCPIII  .  THE  CAPACITY  OP  ARC  I  BEFORE  ANY  INTERDICTION. 

C  ICC  S  T I  I  I  ....  THE  COST  ON  AN  ARC  BEFORE  ANY  INTERDICTION. 

C  ICLASS  .  CONTROLS  THL  SECURITY  CLASSIFICATION  HEADINGS. 

C  IF  111  UNCLASSIFIED!  121  CONFIDENTIAL!  1)1  SECRET,  I A I 

C  TOP  SECRET,  IS)  CONFIDENTIAL  NO  FORM,  161  SECRET  NO  FORN 

C  IT!  TOP  SECRET  NO  FORN,  181  FOR  UFFICIAL  USE  ONLY. 

C  ICOS  Till  ....  UNINTEROIC  TED  COST  FOR  ARC  I. 

C  INNCST  .  THE  COST  OF  THE  NETMURK  FLOH  BEFORE  A  STRIKE. 

C  IMNCUI  .  THE  COST  OF  THE  NETWORK  FLOW  AFTER  A  STRIKE. 

C  INC  .  TMt  INCREMENT  OF  THE  PROFILE. 

C  INCAPII!  ....  CAPACITY  OF  ARC  I  INTERDICTED. 

C  INCOSI III  ...  INTERDICTED  COST  FOR  ARC  I. 

C  1NCRSF  . USED  IN  THE  ONECUT  SUBROUTINE  IN  THE  CALCULATION 

C  OF  ARC  VALUES. 

C  INFERS  .  A  SWITCH  SET  IN  MNCF  IF  THE  FLOW  PATTERN  IS  INFEASIBLE. 

C  10P  . AN  INPUT  WHICH  CUNTROLS  THE  OUTPUT,  IF  IT  IS  ONE,  ONLY 

C  THE  RESULTS  OF  A  STRIKE  ARE  PRINTED,  IF  SET  EOUAL  10  IwO 

c  the  results  of  the  strike  and  the  entire  network  are 

C  PRINTEO.  IF  EQUAL  TO  THREE  TABULAR  OUTPUT  IS  PRODUCED. 

L  IPRINI  . USLO  IN  THE  OUTPUT  ROUTINE. 

C  IPROF  . A  SWITCH  USED  IN  THE  INNER  PRODUCT  SUBROUTINE.  IF  IT 

L  IS  SET  EOUAL  TO  0  THE  COST  OF  THE  UNIVERSAL  ARC  IS 

C  USEO  IN  THE  TOIAL  NETWORK  COST,  IF  IT  IS  SET  EQUAL 

C  IU  1  THE  SUBROUTINE  INNtR  PRUOUCT  DOES  NUT  USE  THE  COST 

C  OF  THE  UNIVERSAL  ARC  WHEN  CALCULATING  TOTAL  NETWORK 

C  COSTS. 

L  IPKOFL  .  IF  SET  EOUAL  TO  ONE  THE  PROFILE  SUBROUTINE  IS  CALLED 

C  AT  THE  ENO  OF  EACH  DAY.  IF  *  2  AFTER  EACH  CUT  ANu  ALSU 

L  AT  THE  END  Of  EACH  OAT.  IF  •  )  NUT  CALCULATED. 

L  ISAVCP  .  SAVES  (HE  CAPACITY  UF  THL  RETURN  ARC  OF  BST  ARL  • 

t  ISAVE  .  SAVES  BSTARC  SUBSCRIPT  IF  IT  WAS  A  TWU-WAY  ARC. 

C  ISVCPI-*  ....  SAVES  CAPACITIES  IN  THE  ONLCUT  SUBROUTINE. 

L  ISWICH .  SAME  AS  I  TWO. 

I  IIHEIII  .  IIME  TO  REPAIR  ARC  I. 

c  I  I  WO .  IF  BSTARC  WAS  ONE-WAY, IT  IS  EQUAL  TO  ONE.  IF  IT  WAS 

C  TWU-WAY  IT  IS  SEI  EQUAL  10  IWO. 

C  JCAP1  .  SAVES  THE  CAPACIIV  OF  THE  UNI  NT  ERU I LT  E  U  CUSI  PORTION  OF 

C  B  Si  ARC  BEFURE  SIRIKING. 

C  JCAP2  .  SAVES  THE  ADDITIONAL  CAPACITY  AT  INTERDICTtD  COST  OF 

c  B  S I  ARC  BEFORE  STRIKING. 

C  JOBS  .  THF  NUMBER  UF  NETWORKS  TO  BE  PROCESSEO. 

L  KCAPl  .  SAVES  I  HE  CAPACITY  OF  I  Ht  UNI  NT  EK  Ul  C  T  F  U  CUSI  PUHIIDN  t,f 

L  IRSIAC  BEFORE  STRIKING. 

L  KCAP2  .  SAVES  THE  ADUlIIONAL  CAPACITY  A!  INILRUICIEU  CUSI  Of 

C  IBSIAC  BEFORE  SIRIKING. 

C  LOWER  .  LOWER  BOUND  VALUE  FOR  NEIWURK. 

C  LWERIII  .  I  HE  LOWER  OOUNO  UN  ARC  I. 

L  MAN  .  1  HE  MEN  NEEDEO  FOR  HEPAIK. 

C  RAXAAC  .  THE  NUMBER  Of  ARCS  IN  THL  NETWORK. 

c  MAxcsriji  ...  the  cost  at  the  jth  profile  point. 

C  NAXELOIJI  ...  THE  F  LOW  AT  THE  JTH  PROF  HE  PUINT. 

C  MAXI  ........  THE  SUASCRIPT  OF  THE  ARC  WlIH  THE  GREATEST  UPPER  BOUND. 

C  MINCSI  . MINIMUM  COST  FLOW  PAMIRS  BIFURl  I  HE  BEST  ARC  TO  STRIKE  IS 

C  CHUSEN. 

(.  HINCUI  . MINIMUM  COS!  FLOW  PAMIRS  WHILE  IESMNG  FOR  THE  BESl  ARC 

i  TO  STRIKE. 

I  MKCSI  . A  TEMPORARY  EOCATIUN  USED  To  HOLD  THE  COST  OF  A  PROF  1 1  C . 

C  N  .  MAXIMUM  NUMHt A  OF  STRIKE  S/UXY. 

C  SI  .  USED  IN  I  HE  OUTPUT  ROUTINE  10  ELSE  WHICH  WRITE  EG  USE. 

C  NLMI  . US  £0  AS  A  SC  I A  ICO  LIST  IN  I  HE  MNCF  SUBRIIUIINE. 

C  SS I  1 1  . CON  I A I  NS  THE  NODE  NAMES  AEIER  THE  NOOES  HAVE  BEEN 

C  ASSIGNED  THE  I R  UNIQUE  NUMBER. 
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C  NfMi  .  T  Ml  NUMBER  UP  MOM  U  MINIS  DESIRED,  CANNOT  EXCEED  IS. 

C  NR  I S  .  ONI  LESS  THAN  NPPTS.USED  IN  I HE  PROFILE  SUBROUT  INE . 

C  NR  .  THE  NUMBER  OP  SIMMS  PER  OAT  UP  IU  N. 

C  NUMI-A  .  USED  AS  TEMPORARY  STORAGE  IN  I  ME  OUTPUI  SUBROUTINE. 

C  PER  100  .  THE  NUMBER  Of  OATS  I  ME  CAMPAIGN  IS  TO  LAST. 

C  PKIII  .  THE  PROBAB I L I  I T  UP  KILL  ION  ARC  I  ASSOCIATEO  MllH  THE 

C  CHOSEN  ATTACK  UNI T. 

C  PPIS  . THE  SANE  AS  NPPTS  BUT  IN  MEAL  MODE. 

C  REPAIRIII  ...  REPAIR  COST  FOR  ARC  I. 

C  RE  PR  I  I  I  .  THE  REPAIR  COST  OP  AN  ARC  AMLR  INTEROICIION  IS  OON  t  . 

C  RESTOR  . .  THE  DAY  ON  WHICH  ARC  I  MILL  BE  REPAIREO. 

C  SAVCAP  . A  LOCATION  USED  10  STORE  THE  CAPACITY  OP  ARC  MARI  ANO 

C  BSTARC . 

C  SAVCST  .  A  LOCATION  USED  TO  STORE  THE  COST  UP  ARC  MAXI  AND 

C  BSTARC. 

C  TITLEIAl  ....  THE  ARRAY  CONTAINING  THE  TITLE. 

C  TUI  I  I  .  II  AT  INPUT  TINE  CONTAINS  THE  LOCATION  OF  I  HE  SUCCLSSOM 

C  NODE  IN  GEOREP  OR  UTM. 

C  it  APTER  NODES  ARE  NUMBERED  IT  CONIAINS  I  HE  UNIQUE 

C  NODE  NUMBER. 

C  TOMAY  .  IP  SET  EQUAL  TO  l>  ALL  ARCS  MILL  61  CONSIDERED  UNE-MAY. 

C  IP  SET  EQUAL  TO  i,  ARCS  WITH  RETURN  ARCS  MILL  BE 

C  CONSIDERED  TWO-WAY,  ALL  OTHERS  MILL  BE  ONE-WAY. 

C  UPPER  I  1 1  ....  THE  UPPLR  BUUNO  OR  CAPACITY  OF  ARC  I. 

.  XLST  .  AN  ARRAY  USEO  AT  INPUT  TIME  CONTAINING  SWITCH  HOMOS. 

C 

c ••••••« •♦••♦•**•«••»♦••••**•*•••**•••«*•♦••••*•♦•••***••*•••♦•»* 

c 

C 


CONNUN  /BLKO/  Ml NC S I t Ml NCO I >N  ,  I  NPE AS 

COMMON  /BLK  1  /  IARI.  PERIliO.  DAY,  MAKARC,  BSTARC,  ICLASS,  IUP,  l  INI  ,S, 

C  BSTARC. 

IMAXNUE  i  JCAPI  ,  JCAP/.Kt API  ,KL AP2 , NR , I RNC ST  ,  I MNCUT  ,  I OMAY 
COMMON  /BLK//ARCNMIBOOI  .  FROM!  800  I  .  TUI  800  I  ,  I  COST  (BOO  I  , 

1 1C  API  dOOl  .  INCAPI  BOO  I  •  L  MLR  I  8001  ,  RE  PR  (  800  I  , 

/REPAIR  I  BOO  I .  I  INI  I  8001  ,FlUw(  8001  i  P  1 1  800  I  ,  NL  I  800  I  .UPPER  I  BUD  I , 
HSlBOOl.NMBUOI.Rt  S  TOR  I  BOO  I  ,  I  CCPI  BOD)  ,  ICCS II  BOO  I  ,  PR  I  d  JO  I 
CUNMUN  /BLK*/  MAXPLUt/SI  .MARCS!  I/Sl  .NPPI  S.  JllMPT  JO)  ,1  IILE  liul 
COMMON  /BLK  S/  I  T  MU.  I6STAC  .  I  SAVCP.  ISA  VC  I  ,  ISmICH.  IPHTiPL.MAN 
INTtGCR  PERIOO.UAY.pl OR. PROM.TU.PI. OP PER. BSTARCiS. REST OR. RT PAIR, 

I R  EPK . ERROR . TOmAY . ARCNM . UC  OS  T 
C 
C 
C 
C 

C  BR INC  IN  DA  I A  NOR. 

C 

CALL  INPUT 
C 

C  NON  NUMBER  TMt  NODES. 

C 

call  number 

c 

C  HERE  Ml  RIPE  UU I  THl  NO JE  NUMBER  SCRATCH  USI  ANO  IHt  PLUM 

C  ARRAY. 

C 

DO  *0  M  »  l.MAXARC 
UPPER  I  Ml  *  0 
ISIMI  «  0 
RE  S TOR  I M I  «  0 
NKMI  *  0 
PLOMIMI  •  o 
*0  CONTINUE 

OL  SCR  THE  IHt  CAMPAIGN  HtKE. 

CALL  OUTPUT (II 
C 

C  SET  UP  A  PLOW  FOR  OUTPCl  COSI. 

C 

CAll  MNCP I  MAXNOt  .  MA  XARC  .FR  OM.  It),  ICUS  T  ,  I L  AP ,  LMI  R  ,  F  l  Ow ,  P  I  ,  NL  ,  INF  E  AS  I 
IF  ( INFERS  .IQ.  II  UU  TO  9*> 

C 

C  WRITE  INITIAL  FLOW  PATTERN  HERE. 

C 

r<  A  Y  *  0 

.ALL  OU I PUI I  6 1 
C 

C  HERE  HE  GET  I  HI  COST  OP  FLOW  IN  IME  NL  T  WORK  (OR  OUlPUT. 

C 


r»  o  n  o  o  r*  non  o  o  o  normrcrnon  o  n  n  r>  c>  c>  c*  r*  n  r*  n  r  n  r  no 
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CALL  INPHDIFIOH'ICOST.MAXARCiIMNCSI , I  ART,  1 1 

call  c/utputi  m 

DAT  -  I 
NR  •  0 
MAN  •  0 

C 

C  START  Next  STHlXE 

c 

Td  NR  ■  NA  ♦  1 
INFCAS  ■  0 

IF  ASKED  FOR  n  PLOT  A  PROF  ILL. 

IF  I  IPAOFL  .£0.  21  CALL  PROFIEU) 

GO  10  ONECUl  10  FINO  IMt  ARC  tO  STRUE. 
CALL  ONECUT 

IFI INFCAS  .£0.  II  GO  TO  95 


BS I  ARC  IS  THE  ARC  10  STRIKE. 


SAVt  TM£  CAPACITY  FOR  OUTPUT. 

KCAP1  •  ICAPI IBS! AC  I 
KCAP2  -  ICAPI IBSTAC  »  II 
JCAP1  -  ICAPIBSTARC I 
JCAP2  -  1C AP I BS I ARC  »  II 
ISAVL  -  B  S  T  ARC 
K  ■  IBSTAC 

WL  INCRtASt  THE  NUMBER  OF  M£n  NEE0E0  FOR  REPAIR  HERE. 

NAN  ■  MAN  ♦  REPAIRIKI 
T5  CONTINUE 

RESTORER  I  ■  OAV  *  ITNEUI 
RESTORIK*!  I  «  RE  S  TOR  I X I 
ISIKI  •  I  SI  R I  *  1 
ISIRMI  ■  I  SIR* 


ACTUAL  INTERDICT ION  OONE  HERE. 


UNINTEROICTIO  COST  ARC. 

J  -  I  SIR  •  ♦  1 
ICAPI R I  •  INCAPIX I 

INCAP  I R I  •  II.  -  PRIKII**J*FLOATIICCPIRI  I  *  .5 
REPAIRIRI  »  FLOATI RE  PR  1X11*11.  -  PX IX 1 1 ••  I  J- 1 1 *PXI  X »  *.5 

INTERDICTS!)  COST  ARC. 

ICAPIR*!)  -  INCAPIX*!) 

INCAPIX*!  )  -  II.  -111.  -  PXIXII)**JI AFLOAT  I ICCPIK*  1 ) )  *  .5 

IF  NOT  T MO* HAY  SHIP  AROUNO. 

IFI  ISMTCH  .EO.  II  CO  TO  90 

PRINT  THE  RESULTS  OP  STNIR1NC  THE  RETURN  ARC. 

CALL  OUTPUT  I  91 

NON  CO  STRIKE  ISTARC. 

X  ■  I  SAVE 
ISMTCH  ■  1 
GO  TO  TS 
90  CONTINUE 
C 

C  MAI  TE  OUT  THE  RESULTS  OP  THE  STRIKE  HERE. 

C 

CAU  OUTPUT  1 21 

C  IP  THE  PLOM  IS  LESS  THAN  THE  OE RAMOS  SAY  SO. 

C 
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95  CONI INUE 

in  inf i  as  .to.  1 1  oo  10  ino 
ififluwiii  .ci.  lweri  1 1 1  co  ro  190 

If  IH1)«I  II  .£C.  0  I  CO  IU  200 
IFIHSIAHL  .lb.  II  GO  10  100 

C 

C  CHICK  10  SO  IF  OAT  IS  00N1 

c 

IFINA  .it).  N  I  GO  10  100 
GO  IU  70 
C 

C  O-fCK  10  SO  IF  CAMPAIGN  IS  DONE 

C 

100  IFluAT  ,t0.  PL  R 100  I  GO  TO  210 
C 

C  IF  ASKtU  FOR,  PLOT  A  PROFIU. 

C 

IF)  IPROFL  .Nf .  II  GO  1C  109 
IFI  I0P  .CO.  21  UNC  •  0 
IFI  I0P  ,NC  •  21  l  INt  -98 

109  IIIIPMOFL  .it.  21  CAll  PROFLE  1 1 1 

C 

C  OPOAI t  FOR  NtXf  OAT 

C 

1 10  OAT  •  DAT  *  I 
NR«  0 

C 

C  II  AN  ARC  IS  TO  8t  MtPAIRtO,  RCSIOKC  AIL  2LR0  SIRlRC  VALUtS. 

C 

00  160  N  •  S iMAXARC  , 2 
IMRtSlORINI  .G1.  OAT  I  GO  TO  160 
IM  I  SI  Ml  .(0.  Ul  CO  TO  160 

C 

C  1NIIIALI7E  I  ML  UNI NTE ROIC  1 1 D  COSl  ARC. 

C 

ICAPIMI  •  ICCPIMI 

INCAPIMI  •  II.  -  PKIMI)*F10AIIICCP(MII  *  ,S 
RLPAIRINI  •  PKIMI  •  FLUAlIRtPNINI  I  ♦  .S 
ISIMI  •  0 
AESIOKIMI  •  0 
160  LClNl INUt 
MM  .  S  ♦  1 

00  I/O  M  •  MM, MAXAKC  ,2 

C 

C  IN  I  I  I AL  I  2  E  THt  INTtHOICIIO  COST  ARC. 

C 

lFIRtSIUHIMI  .Cl.  DAT  I  CO  TO  l  TO 
IFI  ISIMI  .IQ.  01  GO  TO  1  TO 
ICAPIMI  •  0 

INCAPIMI  •  II.  -II.  -PRIM-111  I«FIUATI  ICCP  I M I  I  ♦  .S 
ISIMI  .  0 
SCSI UR (Ml  *  0 
1  TO  CONI INUt 
C 

C  St  I  UP  FLOW  FOR  THt  NEXT  OAT. 

C 

CALL  MNCFIMAXNOt  ,  MAXARC  ,FRUM,  10,  ICOS  1 , 1C  AP,  LWtR  ,  F  LOW  ,  P  I  ,NL,  INF  t  AS  I 
IFI  INFEAS  .EU.  II  CO  TO  9S 
C 

C  GET  I  IS  COSl. 

C 

CALL  INPROIFLOW.I COST  , MA XAKC , IMNC ST  ,  I AR T  ,  1 1 
C 

C  PRINT  NETWORK  FLOW  PAT  URN  AND  IIS  COSl. 

C 

CALL  OUTPUT  I  61 
CALL  OUIPUTITI 
GO  TO  TO 

190  CALL  OUTPUT! II 
GO  TO  100 
200  CALL  OUTPUT  1 61 
GU  TO  100 

210  IF  HOP  .NE .  2)  LINE  •  6« 

IF  IIUP  .Nt.  21  CALL  OUTPUT  161 
CALL  OUTPUT  I Tl 
C 

C  IF  ASKEO  FOR,  PLOT  A  PROFILE. 

C 

IFI  IPROFL  .It.  21  CALL  PROFLilll 

CALL  tXIT 

ENO 


I  HP  HI, 


'.UnHllullM  INPHOIH  Ull.LCU  I  ,MAXARC,HIN(.  I  i  I  AM  I  ,  IPRUM 
I  l  OKI  l  »  .C.USII  l  I 
IMK.t*  KOI)  .  CUSI 


*r<n  -  o 
i  ip  i  <  i  •  i am  r  «  i 

u  i Pxtif  i.  icuai  ru  u  i ml n  ust  ihi  cusi>  ns  tm  univuis/u.  ahc. 

Ill  IPMuf  .fw.  U  I  I  IHI  •  1 
00  10  I  *  IIM[,*A«AKC 
•  INCI  -  minCI  ♦  H0HI1  »»t<)Sl  I  !  I 
0  (  I  OH  N 

£  So 
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INPul 


SUIlRUUflNl  INPul 

0  I  Ml  Si  I  UN  IlHPI/Ol  ,AISM/I 

lOMNUN  /dll  0/  MINI  if  ,  Ml  NCUI  t  N  •  I  III  l  AS 

llJMMUN  /dim/  I  AX  I  ,Pl  MlUU.UAY  .NAAAAI  ,  bil  AMI  ,  I  Cl  Aii  ,  I IJP  ,  l  I  Nt  ,  i  , 

I MAXNOC  , Jl AP I  ,JIAP/,XCAP1  , X(  AP  /  ,  NX  ,  I  MNl.  i  f  .  I MNt UT . I UfcAY 
lOMMUN  /dLM//AHlMN(  800)  ,  f  AUM(  800  I  ,  I  01  800  I  •  1C  Oil  I  801  I  t 

I  I  LAP  I  BOO  I  ,  INIAPI  HUU  I  ,l»l  M  I  8001  i«l  PM  1000  I  . 

/H  t  PA  I P  I  HOUI  •  I  IMI  (  BOO  I  ,f  L  OR  4  8001  .PI  I  800  I  .  M  I  BoO  I  ,  OPPl  M  IMOOl  , 
tl  il  8011  ,NNI  800  I  ,Kl  if  OKI  #00 1  t  IllPIBOOl  ,  I  CL  if  I  800  I ,  PI  I  BOO  I 
lUMMUN  /dlAA/  MAXI  IU1/8I  .HAXISI  l/SI  ,NPPI  i,  J 1 1  NPI  10  1  ,  I  I  til  1/01 
lUNHUN  /  H  l  A  5  /  I  I  ML.  I  Bil  Al  1 1  SAVIP  •  I  5  A  Ml  I  .  I  ihfCH.  IPPIIf  I  .MAN 
IN  IlllM  PIK  lUO.DAT.f  l  OH  ,  I  MUM,  fU.PI  ,UPP  t  M  .  Bi  I  API  .  i  ,  M I  i  I  (IM  ,  Ml  P  AIM  , 
IKlPM.lMKIJM  ,  fUMAT.AMCNM.Ulaif 
INUOlK  IlNP.Xlil 
OAfA  *1 il /• AMI i* , *IUMP'/ 
l 
C 
l 

1*0 

l 

C  MLAO  f  HI  II  III  l  AMU  ANU  SKIP  A  PAOt  . 

I 

Ml  AO  I  b  ,90511  lllll  IMMI, MRU, /Ol 

Ml  ADI  5  ,0001  N.PIMIUO, III Aii, IOP, I  AMI 

1 , IPMOfl ,NPP li, IUWAV, JOBS 
NMIttl  6  ,90AI 
S  •  2  •  I  AH  T  «  I 

l 

C  LOOK  l  ON  f Ml  AMIS  IAAO. 

I 

MIAUI  »  ,90511  ItMPI  JI.JM./Ol 
IMHMPI1I  . In. Al  ill  11  I  00  III  10 
C 

C  II  N0 1  I  Ml  ARCS  IAHO,  IMHOM. 

C 

C  Al L  1X11 

C 

l  Ml  Al)  IN  TMl  AMI  OAIA 

t 

10  I  *  I ♦  1 
l 

C  H  t  AO  AN  UN  I  Nil  MO  1C  It  0  10SI  AMI. 

I 

MIAUI  b  .910IAMCNMH)  .IMUMI  I  1,10111  .IllPIII  .lllill  II.RIPXII  I. 

Ill  Mt  I  II, INI  Mil  I, PAID 
1 1  I AKINM  III  .tU.XLill/ll  00  10  JO 

C 

C  IN  1 1 1  Al  I  !  t  I  Mt  ONI  NIIMOI  C  ft  0  COil  ARC. 

C 

ICAPIII  *  IllPIII 
I  lui  1 1  I  I  •  ICCillll 

INlAPlII  •  II.  -  PAl  I  I  I»FIUAI  I ICCPI  I  I  I  »  .b 
A  IPA  IM  (  |  |  •  PAID  •  f  10AI I  MIPR 1 1  )  I  *  .8 
I  •  I  ♦  l 

l 

1  Mt  AO  I  Ml  COMRt  SP0N01  NO  INIIMOICUO  COSI  ARC. 

C 

MtAOl  b  ,9lOIAHlNM(ll,fHOM(l),IOIII, IllPIII,  ICCillll,  MIPMIII, 

I I  TMl  I  I  I  ,LMIHI  I  I  ,PAI  |  | 

l 

C  INI  I  I  Al  I  It  TMl  INIIMOICUO  10SI  ARC. 

C 

ICAPIII  •  0 
ICOSt  III  •  ICCillll 

INCAP  III  •  II.  -II.  -PAII-llll*FI.OATIICCP  MM  •  .i 

c 

GO  IU  10 
JO  CONI INUl 

NAXAHC  >1-1 
RETURN 
C 
C 
C 

900  FORMA?!  lAIJA.Iltl 
90A  FORMAT! 1H1.10XI 
908  FORMAT  I  JOAN  I 

910  FORMAT  MX,  J  I  AA , JA I , A* , 1 X , M 1* ,AXI ,F 5.A I 
END 


INI  I 
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SURRUUF  Itl  *NC  *  INI. fit  ,,  AMI  >  .  I  ,  J.'.ll'.  I  ,H|  ,ll„H  OK.P  I  ,  NA,  INI  f  AS  I 

i  nncf  is  imc  fuiairson  do i  uf  mui«  au.iiriimn. 

I S r i  M  NOUt  St  ARCS.  I  •  J.odSI  .Hi  ,1  U.HON.M  .NA,  Util  AS 
0  Ml  MS  I  fiN  II  IOC01 . 21  I00GI  ,C(l4  1 1  I'ult  I  .HI  I  I  iUO  1 .1  ill  I ..  I  i  I ,  FI  Utt  I  JUCO  I 
U I  Ml  NS IUN  Ml  Ilium  ,N*I  lOolll 
l 

c 

l  Off  IN  1 1  llIN  III  CAU1NO  StwUfNCl 

c 

c 

C  NANI  USl 

c 

l  NOOIS  NUMNIN  Uf  NOC/fS 

l  Dll  NUNHIR  Of  ARCS 

C  I  list  n  FROM  XM  t  C  I  NM  I  NO  <  NliOt  s 

C  2  list  Of  IU  XI  NO  I  No  <  NOlll  S 

C  CUSI  UN  I  I  COST  Of  flow  CM  ARCS 

C  HI  uPPf*  IIOONUS  FOR  ARC  1 

l  IU  lUtfiK  8UUN0S  fU«  ARC j 

C  flJW  ANOtjN  I  Of  HUM  IN  ARCS 

C  PI  NUDl  PRICfS 

C  NA  SCRAICM  LISI  fuA  NOOt  lABUlNC 

C  INffAS  f  l  AO  OfNIlFING  tHf  IfMOUluV  ill  OufPllf 

c 

C  BII.IN 

IN  If  Of  A  A.AA,N.$*C,SN«,llf  t  ,  IN*  .(  ,  AO*  ,C O*  ,N  l  ,  N2  ,  INC,  l  ABll 
INf  •  2IAfNAI«.Nf 
AO*  •  0 

C  IUO*  PUR  AN  dUI  Of  KIITCR  ARC 
C  OU  90  AA>I,  ARCS 

AA  •  | 

100  N I  •  1 1  A  A  I 

N  2  •  J I  A  A  I 

C  *  LUSIIAAI  'MINI  I  -  PI  IN2I 

NO  II  IHOMlAAI.il  .101  AAI  .UR.  1C.  I  T.O.  ANfi.  f  l  Ow  I  AA  I  .  1 1  .HI  I  AAI  I  Ir.010 
If  IflUMlAAl  .OI.HlI  AAI  .L'R.  If  .01 .0.  AND.  II  OR  I AAI  .01  ,1  UIAAIII  GOIf) 

C  VO  CONIINUL 
9U  AA  >  AA  *  l 

If lAA.lf .ARC  SI  00  TO  100 
C  NO  UUI  Of  AHUM  ARCS  Cf f  I 
INffAS  •  0 
Rf I URN 

C  UUI  Of  KILIM  ARC  FUUNU 
SO  SRC  •  21  AA  f 
SN*  •  | I AAI 
CABU  «  »AA 
00  IU  200 
BO  SRC  ■  1 1  AAI 
SN*  •  >  I A  A I 
lAtlfl  •  -AA 

C  SAVl  LABflS  If  CASI  lIRfKAIION  WAS  INCRIASINC  NO  Of  PRICfS  UN  FHIS  ARC 

200  If  I  AA.fO.AOK.ANO.NA) SRC  I .Nt.OI  GC  FU  2US 
00  201  N  •  I.NUOtS 

NAINI  •  0 

201  CONIINUf 
AO*  *AA 

20)  CO*  •  C 

NAtSMCI  •  l  ABU 
C  L  ABU 

210  l  ABU  •  0 

JO  2S0  A  "I. ARC  S 
N!  •  1 1 A I 

IF  INI.LF.OI  r.U  10  2S0 
N2  •  21  A) 

IF  INAINl).tG.0.ANU.NAIN2l .10.01  OU  FU  2S0 
If  I NAI Nl I .Nt .0 .ANU. NA I N2 1  .Nf  . 0 1  GO  FU  2N) 

C  •  COS  FI  A I  *  FUND  -  RIIN2I 
If  INAINll.tG.QI  GO  FO  220 

IFIFIONIAI.GI.HIIAI.UR.I FLOW IAI.Gf.iniAI.ANU.C.GT.011  OU  FU  2NS 
NAIN2I  ■  A 
GU  TO  2A0 

22U  If  I  FLONl  A  I.Lt.LOIAl.OM.I  F  LUWI  Al. If. HUAI.  AND. C. IT  .Oil  GO  IU  2NS 

NAINI)  •  -A 
2N0  LABIL  ■  I 

C  NUOf  lABilfUt  IfSI  FOR  URIAH  THAU 
IF  INAISNKl. Nt.OI  GO  TO  2B0 
2N)  II At  •  -Nl 
2)0  CONTINUf 
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C  GO  BACK  AMO  DU  MORE  LABELING  If  SUNi  NOOf  MAS  LABELED  ON  CAST  PASS 
If  (LABEL. NE. 01  CO  TO  210 
C  RESTORE  POSITIVE  SIGNS  TO  FIRST  NOOf  LIST 
260  00  2 TO  A  •  I, ARCS 
1 1  A I  ■  I AB  SI  1 1  A  1 1 
2 TO  CONTINUE 

C  IF  NO  LABELING  OONE  UN  LAST  PASS.  GO  TO  INCEASE  PIE 
If  ILABll.fQ.OI  GO  TO  600 
C  BREAK! WU .  fINO  THE  INCREMENT 
100  INC  •  INF 

C  FULLUm  PATH  BACK  FROM  SOURCE 
N  •  SNC 
110  A  •  IABSINAINII 

If  INAINI.LT. 01  GO  TO  IIS 
N2  •  1 1  A  I 

C  •  COS  1 1 A I  -  MINI  ♦  PIIN2I 
If  1C  .Gt  .01  INC  •  HINDI  INC, LOIAI-FLOMIAI  I 

If  IC.LE.OI  INC  •  NINOI  INC.HIIAI-f  LOMIA)  I 

GU  TO  140 
IIS  N2  •  J I A  I 

C  •  COST  I A I  ♦  PIINI  -  PIIN2I 
320  If  1C. LI. 01  INC  •  NINOI  INC.fLOMl  AI-HI  I  Al  » 

If  IC.GE.OI  INC  •  NINOI  INC.fLOMl  AI-LOIAI  I 

340  N  •  N2 

If  IN.NE.SRCI  GO  TO  310 
C  INCREMENT  ARCS 
3S0  A  •  IABSINAINII 

IF  INAINI.LT.OI  GO  TO  360 
FLOMlAI  •  fLOMIAI  ♦  INC 
N  •  1 1  A  I 
GU  10  JTO 

360  FLOMlAI  •  FLOMlAI  -  INC 
N  •  J I A  I 

310  If  I  N.NE.SRC  I  GU  IU  3S0 
C  fLOM  INCREMENTED.  RETURN  TO  KILTER  IEST 
NAINI  •  0 
GU  IU  100 
C  CHANGE  PI 
400  OEL  •  INF 
C  FINO  INCREMENT 

00  420  A-l.ARCS 
Nl  -  1 1 Al 
N2  •  J I  A I 

If  INAIN1I.EC.0  .ANU.  NAIN2I.E0.0I  GO  TO  420 
IF  INAINll.Nt. 0  .ANO.  NAtN2l.NE.OI  GO  TO  420 
C  •  COSTIAI  .  PI  I  Nl  I  -  PIIN2I 

IF  INAIN2I.E0.0  .ANU.  FLOMI A I . L I . HI  I Al  I  OEL-MI NOI OEL . C I 
IF  INAIN2I.NE. 0  .ANU.  F  LOMl  A I  .GT  .  LOI  Al  I  OEI.»M|  NO  I  Otl  . -C  I 
420  CONI INUL 

IF  IOLL.NL.INFI  GO  10  430 

IF  IFLOMI  AAI.CQ.LOIAAI  .OR.  F  L  OM I  AA I  .  E  0.  H 1 1  AA  1 1  GO  TU  42S 
l  INFEASIBLE  SOLUTION 
IMCAS  •  I 
RL  TURN 

C  INCREASE  PI 

42 S  DEL  •  IABSICOKI 
430  00  4S0  N  *1. NODES 

IE  INAINI.EO.OI  PIINI  •  PIINI  *  OEL 
4S0  CUNI INUE 

C  GU  BACK  TO  KILTER  TEST 
GU  TU  100 
END 
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N  001  NO 

fUM  I  ION  NOOENUI  111,  SN t 
U I  ME  NS  ICN  NNI  1 1 

COMMON  /mul/  llll  ,Pt«IUU, OAT, MAXARC, BSIARC,  ICIASS,  loP.L  INI  ,S, 
IMAIN01  ,  JLAM  l  ,  JCAP2.XC API  ,KC AP/  ,NR,  IMNIST  ,  I MNt Ul , E UNAt 

t 

C 

C 

NOblUi;  *  MX  X NOt  ♦  I 

IHmAxNOC  .EO.  01  KElUAN 
00  10  I  -  I.  MAXNOE 
■  FI  III  .to.  NNI III  NOOf NO*  I 
10  CONTINUE 
RE I  WIN 
fcNO 


NUMBE A 


SU8RUUIINE  NUMBER 

COMMON  /BLAU/  MINCSI iMINCUf«N,INHAS 

COMMON  /BlMl/  |ARI,PERIOO  ,OAY  ,  MAX  ARC  i  AST  AKC  ■  ILL  ASS,  101*,  L  INI  ,St 
1MAAN0E  «  JC  AP  l  ,  JC  AP/ ,  AC  AP  l  , XL  AP/  ,  NR ,  I MNL  SI  ,  I MNCU  I  , I OM AV 
COMMON  /BLR2/ARCNMI8GOI  ,  TRUMI  8001 , 101  BOO  I  ,  I CUS!  I  #00  I , 

I ICAPI800I,  INLAPIBOOI  .EMER I  BOO  I  ,R£PR 1 8001 , 

/RE PA  IP  I  8001  ,  I  THE  I  8001  , F  LGWI  BO 01  ,P  1 1  8001  ,  NL  I  BOO  »  ,UPPE  R  1800  I , 
IISIBOOl.NNI  Bool »Rl$  IGA  I 8001  ,  ICCP 1 8031 , ICC S I  I  BOO  I  , PA  I BOo I 
COMMON  /BLAB/  1 1  MO,  I  BSF AC , I SAVCP,  I SAVC I , I  SMI  CM, I PROFt .MAN 
INIECER  PER  100,  OAY.r  LON,  F  ROM,  TO.PI,  UPPER,  AST  ARC,  S.KtSTOR,  ME  PAIR, 
IREPR  ,  ERROR  •  I  ON  AY  ,  AKCNM,  OCOS  T 


C 

L 

C  NUMBER  I  HE  NOUES 

C 

MAXNOE  >  0 
00  10  t-l, MAXARC 
NNI MAXNOE »l I  mf ROM! L  I 
FROMIll  •  NODE  NO  I F  MOM)  L  I  ,  NNI 
MAXNOE  *  MAXOI  MAXNOE  ,  F  ROM|  Lit 
10  CONI  I  NOE 

00  AC  L  *  l.  MAXARC 
NNIMAXNOEMI  *  101 L  I 
lull  I  *  NODE NO I IUIL I  ,NNI 
MAXNOE  «  MAXOI  MAXNOE  •  FOIL  I  I 
AO  CON? I NOE 
C 

C  CHECK  FOR  OE  AO  NODES 

C 

00  A 6  I  ■  | (MAXNOE 
00  A/  t  •  l, MAXARC 

if  ifrumili  ,eu.  it  oo  to  aa 

A/  CONIINOt 

MRIIEI  6  .9/BI  NNI  I  I 
AA  10NTIN0E 

00  A6  L  •  l, MAXARC 

IF  I  Id  L I  .EO.  I  I  CO  10  A* 

A 6  CONTINUE 

MRIIEI  B  .9101  NNI I  I 
A8  CONTINUE 
RETURN 
C 
C 

c 

925  FORMAT  I IHI .  22HNO  ARC  KCINB  RT  NODE  .  AA/IHII 
9)0  FORMAT  I  INI .  22HNO  ARC  ENUB  RT  NODE  .  AA/IHII 
ENO 


n 
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onecut 


SUBRUUT  INC  OMECUt 

CUM  MOM  /  tUO/  MINCSI  .MINCUf  ,N, INFERS 

COMMON  /BlAl/  IAAI, PERIOD, OAT, FUXARC,BSIARC,IUASS,  IOP,l  INI  ,S, 
IMAXNOE  .  JCAPI  •  JLAPZ.RCAP1  ,KCAPZ,NR,  I  MNC  ST  ,  IMNCUI ,  tUMAY 
COMMON  /BIKZ/AACNMISOOI  ,F  AOMI 100 1  ,  TOI  SOU  I .  ICOST  I  BOO  > , 
IICAPISOOI.INCAPIBOOI.IMERIBOOI.RCPRISOOI, 

ZKEPAiRIBOO)  ,  I  TML  (  <001  .f  LOMl  IOOI.MI  100)  ,NLI»OOI,UPPEA|BOO  I. 

)IS<  <0  01  .NNIAOOI.RISTORI  (001 .ICCRIROOI  1 1  CCS I ((OOt  « PR I <001 
COMMON  /UK*/  MAXFLUIZSI  .MAXCSIIZSI  , NPPT S ,JTEMP( JO  I  ,1  lilt  IZOI 
COMMON  /BLAS/  IlMO,  IBSTAC  ,  I  SAVCP,  I  SAVC  T,  ISmKh,  IPAUFL.MAN 
IN! tots  PLAIOO.OAY.P  lOM.FRUM.TO.PI  , UPPER , <S I  ARC , S  ,RE S TOR, RE  PA  |R . 
1AEPR, ERROR,  IOMAY,  ARCNM.UCOSf 
INTtGCR  DELTA, OECAS* 

C 

C 

C 

C  HIRE  ME  INITIALIZE, 

C  ARCS  FOR  MHICH  THE  NUMBER  OF  STRIKES  IS  IhE  SAME  AS  MHlN  THIS 

C  ROUTINE  IS  INTEREO  MILL  BE  CONSIDERED  ^STRUCK.  ONE  TO  nhICh 

C  AN  ADDITIONAL  STRIKE  HAS  KEEN  AOOED  MILL  RE  CONSIOEREO  STRUCK. 

C 

mincut  •  o 

N INC  SI  •  0 
BSIARC  •  1 
LOMER  •  0 
C 

C  SET  UPPER  BOUNDS  AT  INFINITY. 

C 

DO  10  I  «  S • NA  XARC  ,  i 
10  UPPER  I  I  I  •  .’14  148164  1 
C 

C  START  ElOM  REroRE  ME  SELECT  AN  ARC. 

C 

CALL  MNCI  INAXNOE  ,  MA  KARC  ,F  RUN,  I  U,  I  COS  I  ,  It  AP ,  LmER  i  f  LUM ,  P I  ,  NL  .  INF  LAS  I 

C 

C  IF  INFEASIdLE  F  L  Oh  RETURN. 

C 

IFI  INFERS  .Em.  1  I  RETURN 
C 

C  GET  CUSI  (F  FLJM  IN  IHE  SYSTEM. 

C 

CALL  INPKOtFlUM,ICOST,NAXA*C.NlNCSI  ,  I  AM  1, 01 

C 

C  IF  THE  F  l  DM  IS  LESS  THA'.  THE  STRUCK  CAPACITY  LONER  THE 

C  UPPER  HOUND  *  UTMEKMISL  I  HE  UPPER  BOUND  DUES  NOT  CHANGE. 

C 

DU  TO  I  •  S i RA  XARC  ,  l 

IFIFLOMIII  »  FLUMll*ll  .GT.  INCAPIII  *  INCAPIIMII  GO  TO  10 
INCMSE  «  MA  XOI  OtFlUM  I  I  I  -  i  NC  AP  1 1  I  I  *1  I  COST  1 1  ♦  1 1  -  I CUS  T  1 1  I  I 
UPPER  III  •  I  NCR  SE  •  1TMEIII  ♦  REPAlRIII 
JO  CONTINUE 
J5  CONTINUE 
C 

C  FIND  THE  ARC  Ml  IH  THE  O-lAIESI  UPPER  BOUND,  II  IS  MAXI. 

C 

MAXI  •  I 

00  SO  I  •  S ,MA  XARC , l 

IFluPPERIII  .01.  UPPtRIMAX  II  un  10  40 
GO  TO  SO 
40  MAXI  •  I 
SO  CONI INUL 
C 

C  IF  THE  GREATEST  UPPER  6  .UNO  IS  E DUAL  10  THE  LOmER  ROUND  ME  ARE 

C  UJNE,  IHI  ARC  TO  SIRIRL  IS  BSIARC. 

C 

IF  |  UPPEAIHAXI )  .Nl  .LOMER)  C  TO  60 
C 

C  TEST  lu  SEE  IF  BSIARC  R  «S  TmO-MAV. 

C 

00  S2  I  •S.MAXARC.Z 

I  f  I  FROM  IBM  ARC  I  .EG.  TO!  I  I  .ANO.  TUIBSTARCI  .10.  FROMIII 
I  .AND.  ARCGM I BS T ARC  I  .EU.  A'LNMIIII  GO  10  S4 
5Z  CONTINUE 
SJ  CONTINUE 
C 

C  BSIARC  IS  NOT  TmO-MAY. 

C 


I 
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i  r mo  •  i 
issue  •  bsuhc 
co  ru  56 
c 

C  BSUHC  MAS  TMO-kAV. 

C 

66  I  Tkd  •  2 
I  SSI  AC  •  I 

IMTUMAY  .tU.  II  GO  TO  SS 

C 

C  NUM  STOSS  AMAV  THE  CAPACITY  OF  THE  BEST  ARC  AMO  ITS  BYPASS. 

C 

56  CONTINUE 

ISVCPI  •  1C API  I SSIAC I 
ISVCP2  •  ICAPI ISSTAC  ♦  II 
ISVCPI  •  ICAPI  BS I  ARC  I 
ISVCP6  -  ICAPISSTARC  ♦  II 
ISMICH  ■  1  T MO 
C 

C  SET  THE  UNINTEHOICTEO  AND  INTERDICTED  COST  CAPACITIES  TO 

C  THEIR  STRUCK  VALUES. 

C 

ICAPI ISSTAC  I  •  INCAPI ISSTAC  I 
ICAPISSTARC!  •  INCAPIBSTARCJ 
ICAPI ISSTAC  •  11  -  INCAPI ISSTAC  ♦  II 
ICAPISSTARC  p  II  -  INCAPIBSTARC  ♦  II 

C 

C  SET  UP  THE  PLUM  PATTERN  HUH  THE  BEST  ARC  INTERDICTED. 

C 

CALL  HNCFI MAXNOE  .  MAX  ARC  .FROM.  TO.  ICQST.ICAP.LMER.FLOM.PI.NLi  INFERS  I 
C 

C  IF  INFEASIBLE  FLOM  RETURN. 

C 

IFI  INFERS  .EQ.  I  I  RETURN 

C 

C  HERE  ME  GET  THE  COST  Of  FLOM  IN  THE  NETMOMK  FOR  OUTPUT. 

C 

CALL  INPHOI  FLOM.  ICOST.NAXARC.INNCUT.IARI.il 

C 

C  RESTORE  THE  CAPACITY  OF  THE  BEST  ARC  AND  ITS  BYPASS  SO  THAT  THE 

C  INTERDICTION  CAN  BE  DONE  AT  ONE  PLACE  IN  THE  NAIN  PROGRAM. 

C 

ICAPI  ISSTAC  I  -  ISVCPI 
ICAPI  ISSTAC  »  II  •  ISVCP2 
1C  API  BS  T  ARC  I  -  ISVCPI 
ICAPISSTARC  ♦  II  •  ISVCPA 

C 

L  GET  THE  VALUE  OF  STRIKING  THE  BEST  ARC  THEN  RETURN. 

C 

NINCUT  -  MINCST  ♦  I  LONER  -  REPAIR ISST ARC  I  1/ I TNE  I BSTARCI 
RETURN 

C 

C  SINCE  ME  00  NOT  HAVE  THE  BEST  ARC  TO  STRIKE  ME  SAVE  THE 

C  IM  INTERDICTED  A  NO  INTEROICTEO  COST  CAPACITIES  OF  THE  ARC 

C  Ml IH  THE  GREATEST  UPPER  BOUNO. 

C 

60  SAVCAP  -  ICAPI  MAX  1 1 
ISAVCP  -  ICAPIMAXI  ♦  || 

C 

C  ME  NOM  SEE  MHA I  MOULD  BE  IHE  COST  IF  THE  ARC  MITH  THE  GREATEST 

C  UPPER  BOUNO  IS  SELECTED. 

C 

ICAPINAX  1 1  -  I NC API  NAXI I 
ICAPIMAXI  »  II  >  INCAPINAXi  ♦  II 

C 

C  SET  UP  THE  FLOM  MITH  MAXI  OUT. 

C 

call  mncfimaxnoe.maxarc.from.tu.icost.icap.lmer.flom.pi.nl. INFEASI 

C 

C  IF  INFEASIBLE  FLOM  RETURN. 

C 

IFI  INFERS  .EQ.  1  I  RETURN 

C 

C  NOM  GET  ITS  COST. 

C 

CALL  INPROIFLOM,  ICOST.NAXARC.  NINCUT.  I  ART.  01 
C 

C  CALCULATE  A  NEM  UPPER  BOUNO. 

C 

UPPER  (MAX  1 1  ■(NINCUT  -  Ml  NC  ST  I  •!  TNE  I  MAX  1 1  6  REPAIRINAXII 


r  non 


-33- 


c 

C  If  THt  UPPER  BOUND  IS  STILL  GKEAIIR  THAN  I  Ml  l  OWCM  8UUN0, 

L  NAXI  BECOMES  THt  Bt  ST  AML  ffl  STRIKE . 

C 

If IUPPERINAXI I  .GT.  LOWER)  GO  TO  10 
GO  TO  BO 

10  LOWER  •  UPPER  I  MAX  1 1 
BSIARC  •  MAXI 
BO  CONI  I NUt 

DO  100  I  •  S.MAXAMC.2 
C 

C  If  THt  UPPER  BOUND  ON  AN  ARC  IS  LESS  f  HAN  I  Hf  LOWER  ROUND  OX 

L  NO  CHANGE  ON  THE  UPPER  BOUND. 

C 

I  f  I  UPPER  I  I  I  .It.  LOWER)  GO  TO  100 
C 

L  If  THE  f l uw  is  greaier  than  THE  STRUCK  CAPACITY  also 

C  NO  CHANGE, 

c 

IfltLOWIII  ♦  flOWII*  II  .GT.  INCAPIII  ♦  INCAIMI  ♦  111  GO  TO  100 
C 

C  OTHERWISE  calculate  UtLIA. 

c 

INCRSE  •  MAXOI  Off  LOW!  I  I  -  I  NC  API  1 1  I  *  I  1 1  OS  I  I  I  ♦  II  -  1C0SIIIII 
OECRSE  «  MINOI  ICAPIMAXI  |  -  I  NC  API  HA  X  I  I  ,  f  L  OW  I  MAX  I  ♦  111* 

1 1  ICOSt I NAXI  ♦  ||  -  I C  OS  t I  MAX  III 

OELTA  •  ININCUT  *  INCRSE  -  OECRSE  -  MINCSII  •  ITNEIII  *  REPAIRIIi 

If  OELTA  IS  LESS  THAN  THE  UPPt  R  ROUND.  IT  BECOMES  I  HE  NEW 
UPPLR  BOUND,  OTHIKWISL  NO  CHANGE. 

IflOLLIA  .GE.  UPPER  I  I  I  I  GO  TO  iOu 
UPPER  1 1  I  •  DELIA 
100  CONTINUE 
C 

L  RESIURL  THE  CAPACITY  AND  GO  GE  I  ANOTHER  ARC. 

C 

ICAPI  MAX  1 1  «  SAVCAP 
ICAPINAX 1*11  •  ISAVCP 
GO  TO  ii 
ENU 
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UUIFUT 


iUtWUUflNl  UUIFUIIJI 

COMMON  /BLAU/  N|NC4f*MlNCUT«N,INFtA5 

COMMON  /8LM/  I  AM  T  ,FtMOO,DAV,NAIAAC,aST  AHC,  I  CL  ASS  •  IOF.L  INt  ,5. 
IMAINUt  ,  JLAP  1  «  JLAFJ.XCAFl  ,  KC  AF2  #NA,  I  MNCit  •  INNCUT  , TO*AV 
COMMON  /BLK2/AKCNNI 1001  »F  AQN|  tOO),TOIIOOI,l  COSF  I  800  )  i 
I  I  CAM)  *001  ,  INC  At  <100 1 1 1  Ml*  I  too  I  .AtFAIBOOl  , 
mi>AiM*oui,MNtitooiinouitooi'Mitooi.NMtooi.umHtfOoi, 
mi  *ooi  ,nni  tool,  At  slut  i  toot  ,iccpi  tool ,  iccsr isuoi.fkibuoi 

COMMON  /BLA4/  NAIFIUIJ9I  ,MA(CSTI24I  , NFP T S , JT INF  I  101,1  ItlEUOl 
COMMON  /BLA5/  11*0.  I*51AC,ISAVCP,I5AVLI,I5MIC  .IMVHi  MAN 
INI  toe  A  Ft*  IOO,UA»,FION,FAON,»0,PI  .  UPFt  A  ,  tSI  AAC  ,  5 ,  A 1 5  IOA  ,  Al  PA  I A  , 
lAiPM.tAMOA,  TUMAY.AACNM.UCOST 


IF  IJ  .01.  I  I  CO  TO  ItO 
CALL  PAGtl6,iCLAS5> 

■Aim  6  .tuaii  nut  111.1*1,201 
■  Aim  *  ,90riPlA|00,N 
CALL  PAGl Ik  ,  ICl AS5I 
■Aim  6  ,9041 
LINt  •  0 
At  fUNN 
20  I AAC  •  0 
10  LINt  ■  0 

•  Allll  A  ,9011 
40  I AAC  •  I AAC  *  l 

I F I  I  AAC  .Of.  MAXAAC  I  CO  tO  40 

IF  IFLOMl  IAAC-II  .to.  0  .ANO.  OAT  .Nt.  01  00  TO  40 
IFIFLOMI IAAC-1 I  .to.  0  .ANO.  J  .Nt.  fcl  00  TO  40 
LINt  ■  LINt  *  1 
II  •  FAOMI I AAC I 
IJ  •  rOIIAACI 

■Aim  A  «909IAHCMIIIAAC-II,NNIIII,NNIIJ|,ICAPIIAAC-II,IC0STIIAAC- 
1 1 1  .AtFAIAI  I  AAC*  1 1 ,1  IMt  1 1  A  AC*  1 1  ,LNtAllAAC-l  I  •  F  LOM  1 1 A  AC*  I  I  .AtSTOAl  IA 
/  AC  -  1 1 , 1  S I  I  AAC  -  1 1 , 1 C  AF  1 1  AAC  I  •  ICUST I  I  AAC  I ,  F  LOW  1 1  AAC  I  ,PK  1 1  AAC- 1 1 
IF  1 1  I  Nt  .Nt.  491  GO  TO  40 
CALL  FAGt  I  k  ,  ILL  AS5I 
■Airtl  A  ,9041 
CALL  FAOtlA, 1CLASSI 
OU  10  )0 

40  CALL  AAOt Ik, ICLASSI 
■Aim  k  ,9041 
LINt  ■  0 
At  TUAN 
40  CONTINUE 

IF  IL  INt  »  k  .01.  411  00  10  4)0 
IFIIFAOFL  .tU.  21  CO  10  65 
IF IL INt  .Nt.  01  GO  TO  TO 
65  CONMNUt 

■Airtl  k  ,9001 
LINt  <  LINE  »4 
10  CONTINUE 

IPHINI  •  *51  AAC 

IFIISMTCH  .tO.  II  I  FA  I  Nr  «  It  STAC 

■AIIEI  k  ,9011  OAYiNA,  AACNNI I  FAINT  I ,  NNI NUN1 1 , NNINUM2 I , 
l  AES TO* I  I  A*  INT  I,  HUM),  NOMA,  XAF1 ,  ICAFI I  FAINT  I  .MAN,  INNCUT,  ISIIFAINT  I, 
2FLQMI  II 
LINE  -LINt  ♦  l 
At  TUAN 

110  IFUINE  .Nt.  4tlC0  TO  2)0 
CALL  FACE  Ik,  ICLASSI 
■Airtl  k  ,9041 
LINt  ■  0 

IF  I J  .tO.  51  AETUAN 
2)0  IFILINE  .Nt.  0  I  GO  TO  2«0 
CALL  FACE  Ik, ICLASSI 

210  CO  T0l)20,)20,))0,)40,)20,)50,)k0,)70,420l,j 
120  CONTINUt 

NLMU  •  FAOMI  BSIAAC  I 
NUN 2  -  TOI BSIAAC I 
NUN)  •  JCAF1  »  XAP2 

NONA  >  ICAFIBSTAACI  4  ICAFIBSTAAC  4  l) 

IFI  |0F  .CO.  )l  GO  TO  kO 

IF  I L  INt  4  t  .CE.  401  CO  TO  k)0 

■Airtl  k  ,9101  OAT,  NA ,  AACNMIBST  AAC  I  ,  NNI  NUM1 1  ,NM(  NUN2 1  • 

1 A  C  S  TOA I B  S  T  A  AC  I ,  NUN) ,  X  AF I ,  NUM4 , 1 C  AF  I  B  S  T  AAC  I ,  TAN,  I NNCUT  , 

21 SIBSTAAC)  ,FLOMl  1 1 
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LINE  •  l  IM  ♦  • 

IFIIUR  .id.  II  RETURN 
CALL  RAGE  16, ICLASS) 

WRITEI  6  .4061 
l INI  •  0 

CALL  PACI 16. ICLASSI 
GO  10  40 

no  mini  6  . 4 1 i i  oat 

L  INC  »  l  I NE  *  ■ 

RETURN 

140  WRITE!  *  .4121  OAT 
LINE  ■  l INC  ♦  8 
RETURN 
160  CO  TO  40 

160  WR I  111  6  ,4141  INNCST.FLOWllI 
CALI  RACE  16, ICL A&SI 
MR  I  IE  <  6  ,4041 
LINE  ■  U 
RETURN 
I  TO  CONTINUE 

IFILINEM6  .U.  441  CO  TO  410 
WRITE!  6  ,4161 
IMNRRTS  .01.  141  CO  TO  J4u 
00  160  I  *  I.NRRTS 

aRITEI  6  ,41  T  I  I.NAXFLOI  I  I  .MAXCSTI  II 
180  CONTINUE 
186  CO  TO  60 
190  I  •  0 
400  I  •  I  ♦  I 

IM  I  .GT.  141  CO  TO  184 
Nl  •  I  ♦  14 

IF  INI  .GT.  NRRTSI  CO  TO  410 

WAITE!  6  ,4161  I  , MAXFLUI  1 1  , NAXCSf  I  I  I  .Nl  «NAXFLO(NI  I  ,  NAXC6T  INI  I 
CU  TU  400 

4|0  WRITEI  6  ,4ITI  I.NAXFLUIII.MAXCSTIII 
GO  TO  400 
440  CONTINUt 

NUN  1  «  TOIBSIAACI 
NUN4  <  FRONI BIT ARC  I 
NUN}  •  KCARI  «  KC AR4 

NUN4  •  I  CAR  I  I B  6TAC I  •  ICARIIBSTAC  »  II 
IFIIOR  .Ed.  II  GO  TU  60 
IFILINE  »  8  .GC.  461  CU  TO  410 

WRITEI  6  ,4101  DAT.  NA  ,  ARCNNI I  SST  At  I  ,NNI  NUNl  I  ,  NNI NUN2  I  • 

IRESTOAt IRSTACI  •  NON!  , KC  AR I , NUN4 , ICAR  1 1 64 T  AC  I , NAN,  I NNCUT , 

41 41 164 T AC  I  .FLOW! 1 1 
LINE  •  LINE  ♦  6 
RETURN 

410  LINE  «  48 
CO  TO  180 
C 
C 

c 

400  FORNAIIIH  ,/41X,  1 6H4 TRUCK  STRUCK  ,  l  IX ,  14HT0T  AL  EXPECTED, 4*. 
IISHNAX.  NO.  UNIT4  ARC  , 1 4X ,  4 6MEXRCC TEO  NURSE R  TOTAL/ 

4  1 TX,  46HNANE  ARC  ARC  DAT  ARC  ARC  CAR  AC  ITT ,  4X , 

1  44HFLOM  AT  UNINT.  COST  TOTAL  NETWORK  TINES  EKRECTEO/ 
444H  OAT  NI44ION  OF  ARC  FROM,  4X ,4HTO, 6X , 1 4HI 0  6E  BEFORE, 
61X,  4HAF  TER  , 6X  ,  1 4HB E F ORE  AF  TER.6X  , 46HRCRAIR  FLOW  COST  ARC, 

610H  NETWORK, 

T/S8H  NO.  NO.  STRUCK  NODE  NOOE  RESTUREO  STRIKE  , 
a  44HSTRIKE  STRIKE  STRIKE  COST  RER  OAT  STRUCK, 

4  SH  FLOW I 

401  FURNATI4H  ,  1 1  ,4X ,  1 1 , 1 X  ,11  .1 ,  A4 1 ,6X  ,  1 1 ,4X  ,4 1 IX  ,  1 4  1 , 4X,  41 IX,  1 41 , 
I4X,4I1X,I8I,4X,I1,}X,I8I 

404  FORNATI INI, IOX I 

905  FORNATI  /IH  ,10X| 

906  FORNATI IM  ,  44X.40A4  ///I 

90 T  FORNATI1H  ,14*,  41HCANRAICN  WILL  LAST  FOR  ,|},14H  OATS,  WITH  , 

1 1 1,40N  MISSIONS  RER  DAT.  I 

906  FORNAIIIH  ,44X,iSHCAR.  AT  JNINT . ,  I8X  ,4 1HL0WER  FLOW  AT  TO  SC, 
IUX.41HA00.  CAR.  INT.  FlOW ,4X  .6HSTRI KE / 

41X,  5  IMA  AC  FROM  TO  UNINT.  FLOW  RERAIR  RERAIR, 

11X,  4SHFLOW  UNINT.  NERAIKEO  TINES  AT  INT.  FLOW  AT  INT., 
44X,  TM  SUCCESS/ 

44X,  41MNANE  NOOE  NOOE  COST  COST  COST  TINE, 

61X,  STMBOUNO  COST  UN  OAT  STRUCK  COST  COST  COST, 

TSX, SHRROS. I 

909  F0RNATI2H  .1IA4.4XI,  1 4 . 4K, 14,2 1  4X , 16 1 ,4X , 4 1 14  ,  IX I  ,4X, 14,4* ,  1 4, 

I4X.2I4X.I6I , IX, IS, 6X.F4. 2  I 
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910  FORMA  1 1  IH  ,1X,AH0AT  <11.9*4  M.SSIUN  .11. UH  STKUCF  ARC  .  A  A , 

1 1  AM,  FROM  LOCATION  ,AA,1)H  TO  LOCATION  ,AA,)0H.  IT  MILL  Bf  KLSTOK 
2EO  ON  OAT  ,11, 1AM.  TNI  Ml  S  S I  ON/ 2*  ,  ABHRE  OUCI D  TUI  AL  EXPECTkU  CAP  AC 
JIIT  ON  THE  ARC  FROM  ,19,  AH  UNITS  I.I9.99M  OF  MHICH  MAT  FLOM  AT  TH 
At  UMINTkRDICTEO  UNIT  FLOP  COSTI/IX.AH  TO  ,I9,BH  UNITS  !,I9,96H  Of 
9HHICH  MAT  FLOW  AT  INI  UNI  N  TLROIC  ISO  UNIT  FLOW  COSH.  TOTAL  RtRAIR 
*  COST  FOR  IHf  CAMPAIGN  IS  .  I M, IN, . /2X.92HAN0  THE  IOTAS  EIPICIEO  Nk 
FT WORK  FLUH  COST  PER  DAT  IS  .IR.1RH.  THIS  IS  STRIKE  ,11, JAM  AGAINS 
IT  THIS  ARC.  THROUGHPUT  IS  ,  1 10 ,  IH. ////I 

911  FORMAT I1H  ,  )X, SHAMC  , AX , AHF MOM ,A X ,  2H TO,  )K  , SHC APAC  ITT  ,  2X ,  AHCOST  ,  2K , 

iahflon.tx,  unrepaired  on  oat/i 

912  FORMAT  I  IH  ,///10X,2THFL0H  TOTALLT  STOPPED  ON  DAT  ,  110,///  1 
91)  FORMAT!  IH  , ///10X,  ItHMI  Nl  Ml  UM  THRU-PUT  CANNOT  BE  MET  ON  OAT 

1110, ///I 

919  FORMA T I  IH  ,  ///  )2X,2BHCOST  OF  THk  PRESENT  FLON  IS  ,IB,1H. 
2.9X.19H  THROUGHPUT  IS  ,IR,1H.///I 
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Appendix  C 
SAMPLE  PROBLEM 


For  Illustrative  purposes  a  sample  problem  is  solved  on  the  com¬ 
bined  road-rail  network  of  Fig.  3.  The  sources  are  S  and  S',  and  the 
sinks  are  T  and  T'.  It  is  assumed  that  the  LOC  user  is  attempting  to 
maximize  flow  and  that  the  individual  arc  parameters  are  given  explic¬ 
itly  in  the  input  deck. 


Input  Deck 


SAMPLE  ROAD-RAIL  LOC  NETMORK 


s 

ARCS 

TTSS 

2  1 

TT 

1 

SS 

A 

1000 

l  10  2 

-1000 

1 

0 

0 

0 

.0 

TISS 

TT 

ss 

1000 

•1000 

0 

0 

0 

.0 

SSS 

ss 

s 

1000 

0 

0 

0 

0 

.0 

SSS 

IS 

s 

1000 

0 

0 

0 

0 

.0 

SSS* 

ss 

S' 

1000 

0 

0 

0 

0 

.0 

SSS' 

ss 

S' 

1000 

0 

0 

0 

0 

.0 

T  TT 

T 

TT 

1000 

0 

0 

0 

0 

.0 

1  TT 

T 

TT 

1000 

0 

0 

0 

0 

.0 

T*  TT 

T* 

TT 

1000 

0 

0 

0 

0 

.0 

T*  TT 

T* 

TT 

1000 

0 

0 

0 

0 

.0 

S  A 

S 

A 

a 

10 

A 

2 

0 

.29 

S  A 

S 

A 

6 

12 

A 

2 

0 

.29 

S  A 

A 

S 

a 

10 

A 

2 

0 

.29 

S  A 

A 

S 

* 

12 

A 

2 

0 

.29 

s  c 

S 

c 

A 

9 

2 

2 

0 

.29 

s  c 

S 

c 

2 

b 

2 

2 

0 

.29 

s  c 

C 

s 

A 

9 

2 

2 

0 

.29 

s  c 

c 

s 

2 

A 

2 

2 

0 

.29 

S  0 

s 

0 

5 

7 

9 

2 

0 

.10 

S  0 

s 

0 

A 

10 

9 

2 

0 

.10 

S  0 

0 

s 

S 

7 

9 

2 

0 

.10 

S  0 

0 

s 

A 

10 

9 

2 

0 

.10 

A  • 

A 

B 

T 

9 

7 

l 

0 

.17 

A  S 

A 

B 

A 

6 

7 

1 

0 

.17 

A  B 

B 

A 

T 

9 

7 

1 

0 

.17 

a  e 

B 

A 

A 

b 

7 

1 

0 

.17 

A  T 

A 

T 

9 

7 

9 

1 

0 

.11 

A  T 

A 

T 

9 

9 

9 

1 

0 

.11 

a  r 

T 

A 

9 

7 

9 

1 

0 

.11 

A  T 

T 

A 

9 

9 

9 

1 

0 

.11 

•  c 

B 

C 

a 

7 

A 

2 

0 

.28 

B  C 

B 

c 

9 

9 

A 

2 

0 

.28 

B  C 

C 

8 

B 

7 

A 

2 

0 

.29 

B  C 

c 

B 

9 

9 

A 

2 

0 

.29 

a  t 

6 

T 

9 

7 

A 

2 

0 

.22 

B  T 

B 

T 

9 

9 

b 

2 

0 

.22 

B  T 

T 

B 

9 

7 

A 

2 

0 

.22 

B  T 

I 

B 

9 

9 

A 

2 

0 

.22 

C  D 

C 

0 

T 

6 

9 

1 

0 

.21 

C  0 

c 

0 

b 

A 

9 

1 

0 

.21 

C  0 

c 

c 

7 

6 

9 

1 

0 

.21 

C  0 

0 

c 

b 

B 

9 

I 

0 

.21 

C  T 

c 

T 

a 

10 

J 

1 

0 

.27 

C  T 

c 

T 

7 

12 

1 

1 

0 

.27 

C  T 

T 

c 

8 

10 

1 

1 

0 

.27 

C  T 

f 

c 

7 

12 

1 

1 

0 

.27 

0  T 

0 

T 

A 

6 

A 

1 

0 

•  2A 

0  1 

0 

T 

2 

9 

A 

9 

0 

.2A 

0  T 

1 

0 

A 

6 

A 

) 

0 

.2A 

0  T 

T 

0 

2 

9 

A 

) 

0 

.2  A 

S  S' 

S 

S' 

19 

2 

9 

1 

0 

.18 

S  S' 

S 

S' 

10 

9 

9 

1 

0 

.16 

S  S' 

S' 

s 

19 

2 

9 

1 

0 

.IB 
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s  s* 

S' 

s 

10 

1 

s 

1 

0 

.It 

0  0* 

0 

O' 

II 

2 

s 

1 

0 

.20 

0  0* 

0 

O' 

12 

1 

s 

1 

0 

.20 

0  O' 

O' 

0 

It 

2 

s 

1 

0 

.20 

0  O' 

O' 

D 

11 

1 

5 

1 

0 

.20 

S'C' 

S' 

C' 

12 

1 

1 

2 

0 

.14 

S'C' 

S' 

c 

4 

t 

• 

2 

0 

.19 

S'C 

C' 

S' 

12 

1 

I 

2 

0 

.19 

S'C* 

C' 

S' 

4 

t 

1 

2 

0 

.19 

S'O' 

S' 

O' 

11 

4 

T 

1 

0 

.22 

S'O' 

S' 

O' 

} 

A 

7 

3 

0 

.22 

S'O' 

O' 

S' 

11 

4 

T 

1 

0 

.22 

S'O' 

O' 

S' 

3 

A 

T 

3 

0 

.22 

C'O' 

c 

O' 

4 

4 

6 

2 

0 

.  1  7 

C'O' 

C' 

O' 

3 

7 

6 

2 

0 

.17 

C'O' 

O' 

c 

4 

4 

6 

2 

0 

.17 

C'O' 

O' 

c 

) 

7 

A 

2 

0 

.17 

CM' 

C' 

T* 

10 

1 

T 

2 

0 

.It 

Cl* 

C' 

T' 

2 

S 

7 

2 

0 

.11 

C*  f' 

p 

C' 

10 

1 

1 

2 

0 

.18 

c*  r* 

p 

c 

2 

5 

7 

2 

0 

.It 

O'  P 

O' 

T* 

1) 

) 

A 

2 

0 

.21 

O'  P 

O' 

f  • 

2 

S 

A 

2 

0 

.21 

O'  P 

p 

O' 

1) 

1 

A 

2 

0 

.21 

O'  P 

compute 

p 

O' 

2 

s 

A 

2 

0 

.21 

Output  Listing 


The  output  listing  li  given  on  the  following  five  pages. 


SAMPLE  ROAD-MAIL  LOC  NET MURK 
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w  y* 

••^COQOOO^^M'X^*'-—  ^aw^as,**^^**-***^**-*.-—  NN-*^****-*  —  ^  N 

«uQ  .  »  •  •  i  . . . . .  * . 

^U*0Q3Q003  3  0  3300030  3000  0  000  3  00000000  30  J  3  J  ". 

33i 


►-  0000000000000003  0  DO  3  0  3  0  0  3000000000  300  30 

J  I  «- 

o  m  y> 

-i  o 

u.  ^  u 

•  *~00000^<N4«00<3*»  ?  »  *  *  »#*«***^-»«*-*««4*R**.^^^* 

h  O  ^  O  —  —  -4  M  -*  — 

S  w  O  O 

—  u.  u  •* 

I 


► 

« 

o 


^  • 

<*-  OOOOOOOOOO  JOOOOOOOOOOOOOOOOOOOOOOOOOQOC 

J  z  *- 

—  wO 

•  o 

s;^ 


*0000003000  0  0000  3  300000000000000000  00000  0 


o 

# 


« 

o 


A* 

V*1 

Mt 

s 


o  ► 

wU«A>«000000000000000000000000000000000000000 

a  «  o 
o  3 

R-l| 

oj  O 

or 


*^H-O^*'***'‘»0*0if'000«300003-**0*00000%0  —  CJ^OOO^O 

**»■#-*  rg  ^  fS  -*  -•  -A  -• 

S  —  o 

9  ° 


* 

► 

o 


**000000000000003 3000000000300000030000000 

moi 

SwS 


« 

o 

u. 


* 

-i 

J 

w 

* 

3 

M 

* 

t 

* 

o 


«  00000^»N«i^^^  —  -*  —  -  —  —  —  —  m  *s  ~  i*  ' 

3* 

*  — 

•AJ  *- 
* 


«  O  O  O  O  O  ^  ' 


*  o 

UJ  u 

* 


•  ooooooo* 

NlRi  O  —  — 

>0^0 
•  JO- 
|^U  I 


<p-R-00000« 

*  *  o  o  6  6  © 

•  *«oooooo 

*!u . 


•  N»>  a  •  * . . 

00^/i^hN<^jviQM|4h4jfl».aQUNjNQM/0QUi/lO^Ou^U^Q 

N! 


ir 


l^^(A4MOi4O<fl<MtJfl|RUQ|j»*O^iAiA0  0M^M9ljQUNO»i 


SJ  44  *A  •  N>  N 

aiMMab^44uuOQ*a»»-autotoOQ^ 

4 • 


hfN^tAOOU  JOOOOhiN  r>i- 
•  ••••••••• 


I 


COST  <JP  the  PRtStMT  ALUM  IS  f mm uotcMR u t  IS 


OAT  1  NISSIUN  I  SIMUC*  ARC  S'U*  •  FAUN  COCAIIUN  O'  IU  COCAIIUN  S'.  II  Mill  HI  NtStOMO  (IN  OAT  A.  IHt  MISSION 
MtOUCtD  IOIAC  tAPtCIti)  CAP  AC I  IT  UN  IHt  AAC  FAUN  l|  UNI  I  S  I  II  OF  BHICH  NAT  HUM  Al  IHt  UN  IN  I  tUOIC  I  f  U  UNI  I  FlUB  CUSI  I 
10  10  UNITS  |  4  Ut  BHICH  NAT  HOT  Al  IHt  UNI  Nit  *0 1C  1 1 0  UNII  flU«  CUSH.  IUIAC  MtPAIA  CUSI  tUH  IHt  CAMPAIGN  IS 

AND  IHt  IOIAC  tTPICItO  Nil  BOAT  FCOB  CUSI  PI*  OAT  IS  AW.  THIS  IS  SIAIAt  I  ACAINSI  IHIS  AMC.  IHAUUUHPUI  IS 
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• 

, 

• 

• 

• 

M 

• 

• 

• 

mm 

> 

to 

to 

— 

y 

X 

8 

X 

to 

to 

4 

»- 

X  St 

X 

■S' 

x  y» 

X 

s 

z 

y* 

X 

c  o 

a 

a  a 

o 

c 

3 

a 

a 

to  J 

to 

—  to 

to 

y»  y* 

y* 

A* 

S' 

y 

St 

y* 

S' 

y» 

s*  a  — 

y4 

a 

— 

St  B 

to 

St 

a 

— 

St 

a 

— 

y* 

—  o 

o 

/* 

to  a 

s 

1  to 

3 

St 

to 

3 

y* 

to 

I  -*  X 

a— 

X 

X 

to 

a  -( 

X 

—  X 

X 

to 

y 

to 

X 

•• 

< 

-  o 

to 

— 

3 

— 

3 

4 

3 

to 

to 

to 

— 

to 

48 

to 

*• 

X  -  X 

a 

X 

to 

X 

a 

r  - 

X 

— 

'  X 

to 

X 

a 

X 

to 

4 

a 

t 

*-  —  a. 

x 

a— 

X 

X 

X 

X 

to 

m 

X 

X 

— 

— 

X 

X 

to  ' 

x  m 

X 

X 

X 

X 

2 

X 

X 

X 

X 

X 

X 

a 

X 

3  x 

j 

a 

X 

a 

a 

X 

1  A 

a 

X 

a 

3 

X 

to 

•  to 

a 

A 

to 

a 

• 

to 

3 

to 

a 

• 

to 

a 

• 

#  3 

a 

# 

a 

a 

#  o 

a 

)  *  a 

a 

/ 

O 

a 

X 

X 

mJ 

s 

to 

X 

uu 

to 

2 

X 

f 

to 

X 

X 

•• 

X 

r 

to 

X 

X 

k- 

X 

X 

to  *“ 

to 

to 

to 

to 

to 

<to 

to 

►  — 

to 

►  — 

► 

M 

*•  ■ 

X  O  * 

X 

a 

X 

x  a 

X 

X 

to 

<M 

X 

a 

* 

X 

a  «  3 

• 

a  or 

o 

• 

a  or 

3 

a  x 

(3 

• 

c 

X 

3 

• 

a 

to 

— 

to 

to 

SJ 

i 

— 

to 

to 

48 

4. 

to 

X  ~  * 

X 

X 

to 

X 

X  - 

X 

X 

to 

X 

X 

►- 

a 

c  ■ 

3  X  - 

X 

c 

X 

X 

3  X 

4 

o 

X 

mm 

< 

o 

X 

X 

o 

—  st 

84 

44 

S' 

St 

3  X  3  St 

a 

X 

a 

a  x 

a 

y 

a  x 

3 

*8* 

a 

X 

O 

•y 

o 

-  3  to 

to 

5 

to 

-*  a 

to 

a 

to 

to 

to 

u> 

a 

X 

X 

X 

X 

X 

X 

X 

X 

X 

a 

O  to  X 

to 

3 

i 

X 

to 

X 

to 

o 

X 

to 

a  to 

X 

to 

a 

-  r  - 

z 

X 

to 

S'  mm  X 

M» 

y* 

MB 

X 

to 

St  to 

X 

to 
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0.  ABSTRACT 

A  computer  model  for  developing  and  evalu¬ 
ating  a  targeting  strategy  against  an  op¬ 
posing  force's  lines  of  communication, 
litis  aim  is  to  obtain  the  greatest  reduc¬ 
tion  in  enemy  throughput  and  the  greatest 
time  and  cost  of  repair.  The  network  arcs 
(road,  rail,  or  waterway  segments  or  trans¬ 
shipment  points)  are  characterized  by  be¬ 
ginning  and  ending  nodes,  upper  and  lower 
bounds,  interdicted  and  uninterdicted  unit 
flow  costs,  repair  times  and  costs,  and 
the  probabilities  that  attempted  strikes 
are  successful.  The  model  is  programmed 
in  daily  cycles,  with  the  user  specifying 
number  of  days  and  strikes.  Strikes  are 
targeted  one  by  one.  At  the  end  of  each 
strike,  total  LOC  throughput  and  costs 
are  printed  out;  if  desired,  a  detailed 
status  report  and/or  a  profile  of  total 
flow  versus  user  cost  are  also  output. 

The  FORTRAN  program  is  thoroughly  self- 
documented  . 
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